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provide  further  capability  for  prediction  and  anal- 
y s is  of  hydrofoil  hullborne  seakeeping,  a mathematical  model  and 
computer  program  have  been  developed  to  predict  roll,  sway  and 
yaw  motions  of  hullborne  hydrofoil  ships  in  beam  seas.  Predic- 
tions agree  well  with  towing  tank  data  for  a l:20-scale  model 
of  the  PHM  hydrofoil  crafty.. 


1 


SOMMAIRE 


Afin  d'accroltre  les  possibilities  de  prediction  et 
d'analyse  de  la  tenue  en  mer  des  navires  hydropteres  en  flot- 
taison  sur  leur  coque,  on  a tree  un  modele  mathematique  et  un 
programme  d'ordinateur  pour  prevoir  le  roulis,  le  tangage  et 
l'embardee  des  navires  hydropteres  en  mers  du  travers.  Les 
predictions  s'accordent  bien  avec  les  donnees  obtenues  dans  un 
reservoir  de  remorquage  avec  un  modele  a l'echelle  de  1/20  de 
l'hydroptere  PHM. 


i .. 


iii 


TABLE  OF  CONTENTS 


| 


I 


I 

1 

Notation vii 

1.  Introduction  1 

2.  Mathematical  Model  1 

2.1  Equations  of  Motion  2 

2.2  Hull  Coefficients  3 

2.2.1  Added  Mass  and  Damping 3 

2.2.2  Hydrostatic  Restoring  Coefficient  ...  4 

2.2.3  Exciting  Force  and  Moments  4 

2.3  Foil  Coefficients 5 

2.3.1  Nonzero  Forward  Speed  5 

2.3.2  Zero  Forward  Speed 10 

2.3.3  Strut  Wave-Making  Damping  11 

3.  Computer  Program  12 

4.  Comparison  of  Theory  with  Experiment  12 

5.  Concluding  Remarks  13 

! 

Illustrations  14 

References 24  ! 

Appendix:  Computer  Program  Details  ....  25 

A.  Input  Description  25 

1 

B.  Sample  Input  .....  29 

C.  Sample  Output 31 

D.  Computer  Program  Listing  40 


NOTATION 


A 


jk 
1 j k 
* j k 


'La 


n 


w 

C (k) 
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exciting  force  or  moment 
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rolling  moment  of  inertia 
yawing  moment  of  inertia 
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foil  area 
Sear's  function 
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x - coordinate  of  foil  mid  chord 
time  variable 

wave  horizontal  orbital  velocity 
wave  vertical  orbital  velocity 

component  of  wave  orbital  velocity  perpendicular 
to  the  foil 

coordinate  system  (Fig.  2) 
foil  dihedral  angle 
foil  angle  of  attack 
flap  deflection 
rudder  deflection 

two-dimensional  section  potential 
wave  amplitude 

flap  control  system  damping  ratio 
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wave  elevation 

sway  displacement 

roll  angle 
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yaw  angle 

variable  of  integration  in  longitudinal  direction 
density  of  water 
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flap  control  system  natural  frequency 
rudder  control  system  natural  frequency 


ix 


1.  INTRODUCTION 


Although  hydrofoil  ships  will  spend  some  of  their  op- 
erational time  in  the  displacement  condition,  little  attention 
has  been  paid  to  the  theoretical  analysis  of  hullborne  hydro- 
foil seakeeping  until  recently.  Indeed,  Ref.  1,  which  treats 
pitch  and  heave  motions  in  head  seas,  appears  to  be  the  first 
published  work  to  address  this  problem.  The  present  report 
describes  a mathematical  model  to  predict  roll,  sway  and  yaw 
motions  of  hullborne  hydrofoil  ships  in  beam  seas;  also  in- 
cluded is  a computer  program  which  applies  to  craft  with  fully 
submerged  foil  systems  arranged  in  either  a canard  or  airplane 
configuration.  This  work  is  thus  a logical  extension  of  Ref.l 
and  together  they  furnish  computerized  procedures  for  predic- 
ting hullborne  hydrofoil  motions  in  the  five  major  degrees  of 
freedom.  Further,  these  programs  are  applicable  to  a wide 
range  of  hull  and  foil  configurations. 

As  in  Ref.  1,  hull  exciting  forces,  added  mass,  and 
damping  are  computed  by  the  usual  means  of  strip  theory,  and 
upon  these  are  superposed  linearized  hydrofoil  terms.  Predic- 
tions agree  well  with  towing  tank  data  for  a l:20-scale  model 
of  the  PHM  hydrofoil  craft.  However,  because  of  the  limited 
scope  of  the  test  results,  one  should  not  base  general  conclu- 
sions on  this  comparison. 


2.  MATHEMATICAL  MODEL 


The  mathematical  model  is  obtained  by  adding  lin- 
earized hydrofoil  terms  to  the  strip  theory  of  Ref.  2.  The 
most  important  assumptions  and  restrictions  are: 

(1)  Ship  response  is  a linear  function  of 
wave  excitation. 

(2)  Ship  length  is  much  greater  than  either 
beam  or  draft. 

(3)  The  hull  does  not  develop  appreciable 
planing  lift. 

(4)  All  viscous  effects  are  negligible  ex- 
cept for  zero  speed  foil  and  strut 
damping . 

(5)  Hull-foil  interaction  is  negligible. 

In  applying  strip  theory  to  a displacement  hull,  (1) 
to  (3)  are  normally  assumed,  but  (4)  is  changed  to  "all  viscous 
effects  other  than  roll  damping  are  negligible",  and  the  effect 
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of  viscosity  on  roll  damping  is  included  at  all  speeds.  For 
hydrofoil  ships,  however,  which  do  not  have  bilge  keels,  hull 
viscous  damping  is  always  negligible  compared  with  foil  and 
strut  damping.  Assumption  (5)  makes  the  problem  theoretically 
tractable  by  permitting  direct  superposition  of  hull  and  foil 
terms  . 

2 . 1 EQUATIONS  OF  MOTION 

Consider  a hydrofoil  ship  whose  length  is  signifi- 
cantly greater  than  either  its  beam  or  draft  and  assume  that 
this  ship  is  travelling  at  constant  speed  U along  a mean  course 
at  right  angles  to  the  direction  of  propagation  of  a train  of 
long-crested  regular  waves  of  frequency  u>  (Fig.  1).  Let  x,  y, 
z be  a right-handed  orthogonal  coordinate  system  fixed  with 
respect  to  the  mean  position  of  the  ship  with  the  origin  at 
the  mean  position  of  the  centre  of  gravity.  The  positive x-axis 
points  forward  in  the  direction  of  motion,  the  positive  y-axis 
to  port,  and  the  positive  z-axis  vertically  upward  (Fig.  2). 
Denote  sway  by  n0,  roll  by  ri  , and  yaw  by  q . 

The  coupled  sway,  roll  and  yaw  equations  are  given 
below,  using  the  same  subscript  convention  as  in  Ref.  2.  Flap 
(3)  and  rudder  (6)  equations  are  also  given,  with  notation 
similar  to  Ref . 1 . 


Sway  : 


Roll : 


Y aw : 


Flap : 
Rudder: 


(A22+m)ri2  + B22n2  + + B^n*  + + A^n,. 

+ B26*6  + C26n6  + A23^  + + C233  + A26*  + B26* 

+ C266 


f2  (1) 


A42^2  + B42n2  + (A44+I4)fi4  + B44n4  + C44n4  + A46\ 
+ B46\  + C46n6  + A43®  + B43*  + C43$  + A46'6‘  + B46* 

+ C466 


F4  (2) 


A62^2  + B62n2  + A64^4  + B64n4  + C64n4  + (A66  + I6^6 

+ “66*6 + C66n6  + A68®  + B66S  + C686  + A66S  + B66S 

+ C6«5 


F,  (3) 


-“S(kiii4  + k»\  + V*’  + 6 + 2?8“66  + “6S 

-",4(kj;ii6  + k*"6  + %V  + 5 + 25«V  + “«6 


0 (4) 


0 (5) 


i, i ill. i hi  iilUU»M*ll>WP<« 


The  A 's,  B 's,  C 's,  and  F ' s are  ascribed  the 
general  form  ^ ^ 

*i.i  • *!j  + *I.i  <6> 

where  A1*  and  A*  denote  contributions  from  the  hull  and  foils 
i J i J 


respectively. 


Interaction  between  hull  and  foils  has  been 


H F 

ignored.  Expressions  for  the  A^  , A^  , etc.  are  given  below. 


2.2  HULL  COEFFICIENTS 


The  strip  theory  used  to  compute  hull  coefficients 
is  obtained  from  Ref.  2.  Since  an  adequate  derivation  is 
given  therein,  only  the  final  results  are  presented  here. 


2.2.1  Added  Mass  and  Damping 


A 

B 

A 

B 

A 

B 

A 

B 

A 

B 

A 


H 

22 

H 

22 

H 

24 

H 

24 

H 

26 

H 

26 

H 

44 

H 

44 

H 

46 

H 

46 

H 

62 


r ac  U A 

[a22d?  “ (i)2b22 

(7) 

£b22d£  + Ub22 

(8) 

. H , JC.  U , A 

A42  £a24d^  u)zb24 

(9) 

b“2  - A24 as  + ua*4 

(10) 

fa22W  - ^[xAb*2  - /b22dd 

L L 

. U2  A 

U)za22 

(ID 

/b22f,df,  + UUAa22  - ^a22d^] 

Li  Li 

+ ^2bA 

a)2  22 

(12) 

r ac  U2  A 

- ^b44 

(13) 

/b44d?  + UaJ4 

Lj 

(14) 

£a24£d£  " ulIxAb24  " £b24d£1 

U2  A 
co2a24 

(15) 

/b24£dE  + b 1 xAa24  - £a24d£> 

L Li 

+ "l»24 
U)2  24 

(16) 

/a22SdS  - ~2lxAb22  + ^b22d?1 

Lj  Li 

(17) 
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/ b22CdC  + U[xAa22  + / a22dU 

L L 

'*245«  - 22|x4b*4  + /b24dC) 

L L 

/b24^d5  + U[xAa24  + /a24dd 
L L 

/a225zd{  - -2XjI>22  + “2  [ xAa 2 2 + ^a22d51 

/b22S2d5  + Ux2a*2  + ;2[xAbJ2  + /b22d5] 

L L 


The  above  integrations  are  over  the  length  of  the 
ship.  In  practice,  the  length  of  ship  is  divided  into  ten  or 
more  sections  and  the  two-dimensional  sectional  added  mass  (a) 
and  wave-making  damping  (b)  computed  for  each  section  using, 
for  example,  the  Frank  close-fit  method.  a22  and  b22  result 


from  sway  motions,  a44  and  t>44  apply  to  roll,  while  a24  and  b24 

are  due  to  cross-coupling  between  sway  and  roll.  Subscript  and 
superscript  A refer  to  the  aftermost  section. 

Note  that  B44  contains  no  hull  viscous  damping  term. 

This  simplification  has  been  made  since  extensive  calculations 
have  shown  that  hull  viscous  damping  is  negligible  in  compari- 
son to  the  viscous  effects  of  the  foils  and  struts. 

2.2.2  Hydrostatic  Restoring  Coefficient 

The  only  hydrostatic  restoring  coefficient  affecting 
lateral  motions  is  C44 , given  by 

C»,  = ft  GM  (23! 

44 

where  A is  displacement  and  GM  the  metacentric  height. 

2.2.3  Exciting  Force  and  Moments 

F»  = pU/(f2  + h2)d£  - i^h2]  (24! 

L 

f"  - Pi(/(f4  + h4)dC  - A41  (25 

L* 


PC(/[C(f 2 + h2)  - i^]d£ 


.U  uAi 

W2} 


where  C is  the  amplitude  of  the  incident  wave  and  the  integra- 
tion is  over  the  length  of  the  hull.  f^  and  h^  are  the  sec- 
tional incident  and  diffraction  forces,  respectively,  given  by 


f,(S)  = g / n exp(k  z'  + ik  y)di  (27) 

z — z w w 

cs 

f^(c)  = g /( yn3  - zn2)exp(ktfz'  + ikwy)d£  (28) 

h2(C)  = a)  / 4>2(in3  - n2)exp(kwz'  + ikwy)d£  (29) 

h^(4)  = a)  / <J>4(in3  - n2)exp(kwz'  + ikwy)d£  (30) 

Cc 


The  integrations  are  performed  over  the  submerged  hull  section. 
n2  and  n3  are  the  y and  z components  of  the  unit  outward  nor- 
mal to  the  hull  at  (£,y,z).  <p ^ and  are  the  two-dimensional 

section  potentials  for  sway  and  roll  oscillations,  respectively. 

k is  the  wave  number,  given  by 
w 


k 

U)2 

as  — 

(31) 

w 

g 

and 

z* 

* Z+hCG 

(32) 

where  hC(J 

is 

the  height  of  the  CG  above  the  waterplane. 

The 

Frank  close-fit  method  may  be  used  to  evaluate  <J>2  and 


2.3  FOIL  COEFFICIENTS 
2.3.1  Nonzero  Forward  Speed 

The  foil  coefficients  are  derived  in  much  the  same 
way  as  in  Ref.  1.  We  begin  by  considering  a foil  of  dihedral 
angle  T and  resolving  its  lift  force  L and  moment  N into  sway, 
roll  and  yaw  components. 


sway  force  = 

-Ls inT 

(33) 

roll  moment  = 

L(ycosT  + zsinT) 

(34) 

yaw  moment  = 

NsinT 

(35) 

Here,  no  distinction  is  made  between  foils  and  struts.  The  fol- 
lowing sign  convention  is  adopted  for  dihedral  and  anhedral  angles: 
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for 

a 

port  dihedral  foil 

of  angle 

r 

DP’ 

Fi  = 

PDP 

for 

a 

starboard  dihedral 

foil  of 

angle 

PDS  ’ 

Fi  = 

~rDS 

for 

a 

port  anhedral  foil 

of  angle 

rAP’ 

Fi  " 

~rAP 

(36) 

for 

a 

starboard  anhedral 

foil  of 

angle 

F AS  ’ 

Fi  = 

r AS 

Denote  by  L,  and  N,  the 
a d 

lift 

and 

moment 

acting 

on  a 

foil  as  a result  of  swaying,  yawing  and  rolling  motions.  Then, 
from  Ref.  1,  equation  (21), 


Ld  - lnc+lc  (37) 

where  the  subscript  NC  denotes  noncirc.ulatory  and  C circulatory. 
In  equations  (23)  and  (24)  of  Ref.  1,  we  substitute 

-n^sinT  + (ycosT  + zsinT)r|^  for  ?. 

H,sinr  for  0 

D 

and  obtain 

lnc  = (f) 2 1 (®n6  - un6)3inr  + fi2sinr 

-(ycosT  + zsinTJn  ] (38) 

Lc  ■ |puscLc,c(k)[((s-f)n6-un6)sinr 

• • L. 

H-r^sinT  - (ycosT  + zsinr)n^  ] - g-^-C(k)yn4  (39) 


where  the  last  term  in  (39)  has  been  obtained  by  intuitive 
analogy  with  the  last  term  of  equation  (24)  in  Ref.  1. 


Similarly,  from  equation  (22)  of  Ref.  1 

Ttpb  c 3 s inT 


“LNCS-LcX- 


16 


-(un6  + fn6) 


(40) 


Consider  now  the  foil  exciting  force  and  moment  and 

denote  by  L„  and  N„  the  lift  and  moment  due  to  wave  .action  on 
W W 

the  foil.  Then  from  equations  (38)  and  (39)  of  Ref.  1, 


lw  - ipuscLase(k)w  + |£c(k)n  (41) 

Nw  - -xLw  (42) 
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where  n is  wave  elevation  at  mid-chord  and  w the  component  of 
wave  orbital  velocity  acting  perpendicular  to  the  foil: 

w = wcosT  + usinT  (4! 


where  w is  the  vertical  component  and  u the  horizontal  compo- 
nent. u is  regarded  as  positive  in  the  direction  of  propaga- 
tion of  the  seaway.  For  beam  waves, 

-k  h ik  y 

w w . , . 

u = we  e (4^ 


-k  h ik  y 
w w 
w = io)e  e 


H = e 


ik  y 
w 


where  k is  wave  number, 
w 

Substitution  of  equations  (37)  to  (46)  into  (33)  to 
(35)  yields  the  foil  coefficients  listed  below.  Summation  is 
over  all  foil  and  strut  elements. 


k\2  = TTPEb(|)2sin2r 
B22  = >£SCLaC(k)sin2r 

= — ffpEb  (y)  2 s inT  (ycosT  + zs  inT  ) 

= _ BCLa^ sinr(ycosr  + zs  ini’ ) 

F i , ^Ct 

c24  = "fpu  Es  TtTc(k)  ysinF 

A^  = A^  = trpEb(-|)2ssin2r 

= pUZsin2r[-irb(|)2  +|sCLaC(k)(s  - |)  ] 


-YPU2zscLac(k)sin2r 

^ ikwy  ^CL 

-ipUZSe  w sinr  [U-g^1  C(k) 


+ Ct  S (k)ue 
La  e 

TipEb(^)  2 (ycosT  + zsinT) 2 


— kyh 


(sinf  + icosf)  ] 
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= ypUESC^C (k)  (ycosT  + zsinT ) 2 
1 , 9cl 

= C(k)y  (ycosT  + zsinT) 

= = -it pEb  (^)  2 ss ini’ (ycosT  + zs ini’ ) 

= pUZsinr(ycosr  + zsinr)[TTb(-^)2 

-iscLaC(k)(S-f)] 

= yPU2ESCLc(C  (k)  sinT  (ycosT  + zsinT) 

i iky  3Cl 

= -^pUESe  W (ycosT  + zsinD[U  gj,  ^(lO 

-k  h 

+ CLaSe(k)e  W w(sinE  + icosT  ) ] 

= •|-PUExSCLc(C  (k)sin2r 

= - ^pUExSC^^:  (k)  sinT  (ycosT  + zsinT) 

1 , 9ct 

= - ■^■pU  2 E xS-g-^i  C ( k ) ysinr 
= 7TPE[s2b(|)2  +^]sin2r 
= pUE[-nb(|)2  +|xSCLaC(k)](s  - J)sin2r 

= -ipU2ExSCLaC(k)sin2r 

i ikwy  9cl 

- -ipUExSe  w sinr[U-^  C(k) 

-kwh 

+ C,  S (k)a>e  (sinr+icosT)] 
Lot  e 


(57) 

(58) 

(59) 

(60) 
(61) 

(62) 

(63) 

(64) 

(65) 

(66) 

(67) 

(68) 

(69) 


The  flap  coefficients  are  obtained  by  using  equations 
(50)  to  (53)  of  Ref.  1 to  evaluate  the  lift  and  moment  due  to 
deflecting  a flap  through  angle  B (Fig.  3).  Resolution  of  this 
force  and  moment  via  equations  (33)  to  (35)  results  in  the  flap 
terms  given  below. 
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IPPT 


A2&  - -2pbFT1(-^)3sinr 

B26  •-5‘,UbF'F(T4-^CLaC(k),il>8l''r 


C2B  " PU2bFcFCLaC(k)e6sinr 


A4g  * 2pbpT1 (-— ) 3 (ypco8r  + zpsinr) 


lPUbFc^T4'  2¥CLaC(k)TLl)(5VCOsr  +zF8lnn  (M) 


C4g  " -pU2bFcpCLaC(k)eg(yFcosr  + Zpsinr) 


A6B  = A2gs  + 2pbp(-jL)‘*(T7  + pTj^Jsinr 


B6B  ■ >VFsl"rI-V+nCLaC<k)Tu’1 


- T -pT,  +iTll)] 


2 ' 1 8 4 2 11- 


C66  = C28X  * 2PU2bF(^)2(T4  + T10)8inr 


where  p is  the  distance  from  the  flap  hinge  line  to  mid-chord 

divided  by  the  semi-chord  (see  Fig.  3).  eQ  is  the  flap  effec- 
ts 

tiveness  parameter.  The  T^’s  are  given  in  Ref.  1.  The  con- 
tribution from  both  port  and  starboard  flaps  has  been  summed  in 
the  above  equations.  y and  zr  apply  to  the  port  flap,  and  8 

r r 

is  positive  for  port  flap  down. 

Consider  now  the  rudder.  Side  force  due  to  rudder 
deflection  may  be  calculated  by  substituting  6 for  <p  and  - — 
for  s in  equations  (23)  and  (24)  of  Ref.  1.  Then  ^ 

Lr  = Trpb(|)2(-|6  - U(S)  +|pUSCLaC(k)(-|(5  - U6)  (79) 


where  L is  rudder  side  force,  assumed  positive  when  acting  in 
R 

the  negative  y-direction  in  keeping  with  our  convention  regard- 
ing dihedral  angles.  Rudder  moment  is  given  by 
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Nr  = TTPb(|)2  s(|  6 + u6)+ ipUSCLc(C(k)x  (|  6 + U6) 

r ^ • c 6 

.irpb£_(U5+_) 


(80) 


A26  _1TpbR  16 


Substitution  of  (79)  and  (80)  into  (33)  to  (35)  yields 
the  rudder  terms. 


(81) 

(82) 

(83) 


‘26 


'26 


-pus^c  + c^coo) 


-|(,°2s8CLaC<k) 


A46  = _A26ZR 


B46  " -B26zR 


'4  6 


06 


66 


“C26zR 


A26SR  + 

-pusRX(”l8R-T1+  claC(k)V 


C66  " C26XR 


(84) 

(85) 

(86) 

(87) 

(88) 

(89) 


2.3.2  Zero  Forward  Speed 

At  zero  forward  speed,  viscous  drag  forces  opposing 
lateral  motions  act  on  the  foils.  By  regarding  the  foils  as 
oscillating  flat  plates  and  equating  the  energy  dissipated  by 
the  non-linear  viscous  effect  during  one  cycle  to  that  dissi- 
pated by  a linear  damping  term,  we  obtain  the  following  viscous 
roll  damping  coefficient: 


B4  4 * ^-P0)f)4E(y2  + z2)  '2  SCnsina 


(90) 


where  n,  is  roll  amplitude  and  C is  the  normal-force  coeff icient 
4 n 

for  a flat  plate  tilted  at  angle  a to  the  flow.  From  Ref.  3, 
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0.0467a  a < 40' 


1.17 


a > 40' 


and  from  geometrical  considerations 


i y / z + tanT 


Similar  equations  may  be  derived  for  the  other  foil 
damping  terms,  but  these  are  much  less  significant  than  the 

viscous  roll  damping  term.  Equations  are  given  below  for  B* 

F 4Z 

and  B . 

6 6 

B22  = In  pa)n2  £SCnsina  (93] 

B^  = -T—  pwn,  ESC  | s | sina  (94] 


where  n and  n are  sway  and  yaw  amplitudes,  and 
Z o 


2.3.3  Strut  Wave-Making  Damping 

A strut  in  or  near  the  free  surface  will  generate 
waves  when  oscillated  laterally.  The  resultant  damping  terms 
due  to  wave-making  affect  roll  significantly  at  low  speeds. 

For  a vertical  strut,  the  sway  wave-making  damping  term  is 

b“2  - f <96; 


U)2b 

where  C,,  is  a function  of  . A curve  obtained  using  the 

W g 

Frank  close-fit  method  is  given  in  Fig.  4. 


Roll  and  yaw  wave-making  damping  terms  are  obtained 
W 

by  multiplying  B_„  by  the  appropriate  foil  coordinates. 


. 


L 


! 


I 
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3.  COMPUTER  PROGRAM 


Based  on  the  foregoing  mathematical  model , a computer 
program  has  been  developed  to  predict  hullborne  hydrofoil  lat- 
eral motions  in  beam  seas.  A program  listing  is  given  in  the 
Appendix,  together  with  detailed  descriptions  of  input  and  out- 
put. Note  that  since  hull  viscous  damping  is  neglected,  this 
program  applies  only  to  the  "foils  down"  case.  A further 
restriction  is  that  the  foil  system  must  be  of  the  fully  sub- 
merged type  and  either  canard  or  airplane  in  configuration, 
i.e.  one  foil  unit  in  an  inverted  T while  the  other  is  either 
an  inverted  ir  or  two  inverted  T's.  Full  details  are  given  in 
the  Appendix. 


>. 

4.  COMPARISON  OF  THEORY  WITH  EXPERIMENT 


The  Davidson  Laboratory  has  recently  measured  wave- 
induced  motions  for  a l:20-scale  model  of  the  220-ton  PHM 
hydrofoil  craft  during  hullborne  operation  in  sea  states  3 and 
5 (Ref.  4).  Representative  wave  height  spectra,  as  measured 
during  the  tests,  are  shown  in  Fig.  5 for  the  full-scale 
craft. 


Unfortunately,  Ref.  4 gives  rather  scanty  lateral 
motion  data  because  of  towing  tank  test  restrictions.  The  only 
useful  frequency  response  measurements  are  for  beam  sea  rolling 
at  zero  speed  (Fig.  6).  Root  mean  square  roll,  yaw  rate,  and 
lateral  acceleration  were  measured  across  the  speed  range  in 
sea  state  3 (Fig.  7),  but  the  rather  academic  nature  of  this 
spectrum  (Fig.  5)  does  not  permit  generalizations  based  on 
these  results,  since  little  or  no  seaway  energy  is  present  in 
the  frequency  range  of  greatest  interest  (.3  to  1.5  rad/sec). 

Fig.  6 shows  generally  satisfactory  agreement  between 
computed  and  measured  beam  sea  roll  response  at  zero  speed.  One 
may  reasonably  conclude  from  this  comparison  that  hove-to  roll- 
ing predictions  should  be  satisfactory. 

Predicted  and  measured  beam  sea  root  mean  square  lat- 
eral motions  are  compared  in  Fig.  7.  Agreement  is  satisfactory 
but,  as  mentioned  above,  because  of  the  peculiar  nature  of  the 
seaway  spectrum,  one  cannot  base  general  conclusions  on  this 
compa  r ison . 
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5.  CONCLUDING  REMARKS 


Although,  as  demonstrated  above,  predictions  agree 
well  with  the  measurements  available,  the  latter  are  not  suf- 
ficiently extensive  to  permit  meaningful  assessment  of  the 
general  reliability  of  predictions.  One  may  reasonably  expect, 
however,  that  hove-to  rolling  predictions  should  be  satisfac- 
tory as  indicated  by  the  agreement  between  limited  experimental 
data  and  predictions. 

Computational  experience  has  shown  that  the  foils  and 
struts  dominate  hullborne  lateral  motions,  even  at  zero  speed, 
and  this  dominance  becomes  more  pronounced  with  increasing 
speed.  Foil  system  damping  completely  swamps  hull  damping,  and 
at  nonzero  speeds  the  dominant  forcing  function  arises  from 
action  of  the  horizontal  component  of  wave  orbital  velocity  on 
the  struts.  Further,  the  control  system  is  effective  in  reduc- 
ing roll  angles,  particularly  for  full-scale  speeds  in  excess 


of  10  knots. 


The  present  work  and  Ref.  1 together  furnish  comput- 
erized procedures  for  predicting  hullborne  hydrofoil  motions 
in  the  five  major  degrees  of  freedom.  However,  the  present 
work  applies  to  beam  seas  and  Ref.  1 to  head  seas.  Work  is  in 
progress  to  synthesize  the  two  and  produce  a computer  program 
which  will  predict  motions  in  five  degrees  of  freedom  at 
arbitrary  headings  to  the  sea. 
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DIRECTION 


DIRECTION  Or 
SHIP  ADVANCE 


SHIP  AND 


E A DIRECTIONS 
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APPENDIX 


COMPUTER  PROGRAM  DETAILS 


The  computer  program  applies  to  a hydrofoil  ship  with 
a fully  submerged  foil  system  of  either  a canard  or  airplane 
configuration.  The  main  foil  is  an  inverted  n (Fig.  8),  while 
the  bow  foil  (canard  configuration)  or  tail  foil  (airplane  con- 
figuration) is  an  inverted  T.  Specification  of  a third  strut 
on  the  main  foil  unit  is  an  optional  input;  another  option  is 
to  split  the  main  foil  into  two  T's.  The  bow  (or  tail)  foil 
also  acts  as  the  ship's  rudder,  and  the  flaps  for  roll  control 
are  on  the  outboard  tips  of  the  main  lifting  foil. 

A.  INPUT  DESCRIPTION 


(a)  ONE  CARD,  FORMAT  (8F10.4) 

U speed  (kt) 

EL  length  between  perpendiculars  (ft) 

HCG  height  of  CG  above  waterplane  (ft) 

XCG  distance  from  CG  to  forward  perpendicular  ( f t ) 

RRG  roll  radius  of  gyration  t EL 

YRG  yaw  radius  of  gyration  * EL 

DISP  displacement  (tons) 

RHO  fluid  density  (slug/ft3) 

(b)  ONE  CARD  (12.2F10.3) 

NFR  number  of  frequencies  at  which  responses  are  to  be 

calculated 

FR1  lowest  frequency  (rad/sec) 

DFR  increment  in  frequency  (rad/sec) 

Notes  (1)  If  computing  motions  in  irregular  seas  with 

U > 0,  set  NFR=18 , FR1-.3,  and  DFR=.l.  If  U-0,  it 
may  be  necessary  to  set  DFR=.05. 

( c ) ONE  CARD  (213) 

NSEA  number  of  sea  states  (maximum  of  10) 

NPOS  number  of  positions  at  which  swaying  motions  in 

irregular  seas  are  to  be  computed  (maximum  of  10) 

Notes  (1)  If  motions  in  irregular  waves  are  not  desired, 

use  a blank  card  for  (c). 

(2)  If  NSEA=0 , ignore  data  cards  (d)  and  (e) . If 
NSEA  > 0,  but  NPOS-O,  ignore  data  card  (e) . 
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NSEA  CARDS  (2F10.4 


HSW(I)  significant  wave  height  (ft) 

TSW(I)  energy-averaged  wave  period  (sec) 

Notes  (1)  Fig.  9,  obtained  using  the  data  of  Ref.  5,  is 

offered  as  a guide  to  the  variation  of  TSEA  with  HSEA. 
Caution  should  be  exercised  in  applying  this  curve, 
however,  since  considerable  variation  of  wave  period 
with  significant  wave  height  is  exhibited  by  natural 
seaways  (see,  for  example.  Fig.  1 in  Ref.  6). 

(e)  NPOS  CARDS  (2F10.4) 

XPOS(I)  x - coordinate  of  position  I (stations  aft  of  FP) 
z - coordinate  of  position  I (ft  above  CG) 

( f ) ONE  CARD  (F10.4) 

GMIN  metacentric  height  (ft) 

Notes  (1)  If  metacentric  height  is  not  specified  on  input, 
(i.e.  GMIN=0),  the  program  will  use  a value  computed 
from  the  offset  data.  This  GM  is  not,  however,  cor- 
rected for  internal  free  surfaces. 


(g)  ONE  CARD  (3F10.4) 

FIAV  expected  roll  amplitude  (deg) 

YAWAV  expected  yaw  amplitude  (deg) 

SWAYAV  expected  sway  amplitude  (ft) 

Notes  ^1)  These  data  are  only  required  for  U*=0.  For  U > 0, 
use  a blank  card. 

(2)  When  computing  motions  in  irregular  seas,  set 
these  inputs  equal  to  1.25  times  the  expected  root 
mean  square  values. 

(h)  ONE  CARD  (12) 

NST  number  of  stations  for  which  offsets  are  input 

Notes  (1)  The  program  assumes  a 20-station  hull  representa- 
tion, with  station  0 at  the  forward  perpendicular  and 
station  20  at  the  transom. 

(2)  The  maximum  value  of  NST  is  25.  However,  since 
the  foil  system  dominates  lateral  response,  in  the 
interest  of  computational  efficiency  it  is  generally 
desirable  to  use  no  more  than  10  stations  to  define 
the  hull.  These  should,  however,  be  equally  spaced 
and  include  the  transom. 

(3)  One  each  of  data  cards  (i),  (j)  and  (k)  is  re- 
quired for  each  of  the  NST  stations. 


(i)  ONE 

CARD 

(F10. 3) 

XA(I) 

station  number 

(j)  ONE 

CARD 

(8F10.4) 

YA ( I , J ) 

J-l. 

8 horizontal  offsets  of  station 

I (ft) 

(k)  ONE 

CARD 

(8F10.4) 

ZA ( I , J ) 

J-l, 

8 vertical  offsets  of  station  I 

(ft) 

Notes 

(1) 

each 

(2) 

Exactly  8 offset  points  must  be  specified 
station . 

The  first  point  is  at  the  intersection  of 

centerline  with  the  station  contour  while  the  eighth 
point  is  at  the  intersection  of  the  load  waterline 
with  the  station  contour  (see  Fig.  10). 

(3)  The  vertical  offsets  are  input  as  heights  above 
hull  baseline  (waterline  zero). 

(4)  The  points  and  the  straight  lines  between  them 
should  provide  a good  geometric  description  of  the 
station  shape. 

(l)  ONE  CARD  (II) 

NSTRUT  number  of  struts  on  main  foil  unit  (2  or  3) 

(m)  NSTRUT  + 3 CARDS  (8F10.4) 

GAM(I)  input  dihedral  angle  (deg) 

SWEEP(I)  quarter-chord  sweep  angle  (deg) 

ALF(I)  angle-of -at tack  relative  to  zero  lift  (deg) 

B ( I ) span  (ft) 

CR(I)  root  chord  (ft) 

CE(I)  tip  chord  (ft) 

TC(1)  thickness/chord  ratio 

Notes  (1)  The  number  system  is  shown  in  Fig.  11. 

No.  1 - main  foil  outboard  tip 
No.  2 - main  foil  inboard  span 
No.  3 - bow  lifting  foil 
No.  4 - main  foil  outboard  strut 
No.  5 - bow  foil  strut 

No.  6 - main  foil  centre  strut  (if  present) 

(2)  The  method  of  inputting  dihedral  angles  is 
shown  in  Fig.  11.  These  angles  are  converted  to  the 
conventional  form  (equation  (36))  internally. 
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(n)  ONE  CARD  (5F10.4) 


X(4)  main  foil  strut  x - coordinate  (ft) 

Y(4)  distance  from  main  foil  strut  tip  to  centre  line  (ft) 

II  ( 4 ) depth  of  main  foil  strut  root  (ft) 

X(5)  bow  foil  strut  x - coordinate  (ft) 

H ( 5 ) bow  foil  strut  root  depth  (ft) 

Notes  (1)  Strut  x - coordinates  are  measured  from  the  quarter- 
chord  line  to  the  CG.  X(4)  is  negative,  X(5)  is  posi- 
tive . 

(2)  Y(4)  and  H(4)  are  shown  in  Fig.  11.  For  the  par- 
ticular case  shown,  Y(4)=B(2).  If  Y(4)  > horizontally 
projected  value  of  B(2),  the  main  foil  is  assumed  to 
be  split . 

(3)  Strut  tips  are  taken  to  be  at  the  intersection 
of  the  struts  with  the  lifting  foils  (i.e.  project 
the  foils  and  struts  through  the  intersection  pods). 

(o)  ONE  CARD  (5F10.4) 


BF  flap  span  (ft) 

PF  distance  from  hinge  line  to  mid-chord  + semi-chord 

EFF  flap  effectiveness 

WF  flap  control  system  natural  frequency  (rad/sec) 

ZETF  flap  control  system  damping  ratio 

Notes  (1)  The  flap  is  assumed  to  extend  to  the  tip  of  foil 

No.  1 (Fig.  11). 

(2)  Flap  effectiveness  is  plotted  against  flap-chord 
ratio  in  Fig.  12.  Note  that  this  plot  is  based  on 
aerodynamic  data  and  that  considerable  doubt  exists 
as  to  whether  flaps  are  as  effective  in  water  as  they 
are  in  air  . 

( p )  ONE  CARD  (2F10.4) 

WR  rudder  control  system  natural  frequency  (rad/sec) 

ZETR  rudder  control  system  damping  ratio 


Notes 


QFDD 

QFD 

QF 


! (1)  It  is  assumed  that  the  bow  foil  is  the  rudder. 

ONE  CARD  (3F10.4) 

roll  acceleration  gain  (sec2) 
roll  velocity  gain  (sec) 
roll  gain 


(r)  ONE  CARD  (3F10.4) 

QRDD  yaw  acceleration  gain  (sec2) 

QRD  yaw  velocity  gain  (sec) 

QR  yaw  gain 
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B,  SAMPLE  INPUT 

A sample  case  of  FHROLL  Input  data  is  given  on  the 
following  page  for  a hypothetical  400-ton  hydrofoil  ship  at 
a speed  of  10  knots.  Note  that  the  hull  is  trimmed  up  1%° 
and  that  offsets  are  given  in  local  section  coordinates,  i.e. 
the  first  point  is  (0,0)  for  all  stations. 


10.0 

18  .3 

3 8 

8.0 

10.0 

12.0 

3.0 

5.0 

10.0 

15.0 

20.0 

10.0 
l 0.0 

10.0 
0.0 


0.0 

.0592 

.9160 

1.377b 

1.6369 

2. 1 8U  1 

2.7553 

3.2105 

0.0 

.7379 

1.0755 

2.213d 

2.9519 

3. 3896 

0.0278 

5.168 

0.0 

2 

0.0 

.8639 

1.7279 

2.6918 

3.0558 

*.3197 

5.  1836 

6.0075 

o.o 

.7970 

1 .5909 

2.392? 

3.1897 

3.987  1 

0.78*5 

5.582 

6. 

2 

0.0 

1.277 

2.550 

3.631 

5.108 

6.385 

7.662 

8.939 

0.0 

.8569 

1.7137 

2.6706 

3.0275 

0.25*3 

5.1*1? 

5.99M 

8.0 

2 

0.0 

1.6291 

3.2583 

0.6870 

6 . 5 1 b6 

8.1*57 

9.77*9 

11.00* 

0.0 

.9172 

1.8300 

2.7516 

3.6688 

* .5860 

5.5033 

6 . * 1 5 

10. 

2 

0.0 

2.2017 

0.0833 

6.7250 

6.9667 

1 1 .208  5 

13. *5 

13.5016 

0.0 

1 .1052 

2.2103 

3.3150 

0 . i*?06 

5.5257 

6.63 

6.831 

12.0 

2 

0.0 

2 • 3233 

0 . b0  6 7 

6.9700 

9.2933 

1 1 . 6 1 6 7 

13.90 

10.10? 

0.0 

1.0627 

2.1650 

3.2080 

0.3307 

5.013* 

6.5 

7.2o7 

10. 

2 

0.0 

e • 3233 

0.6067 

6.9700 

9.2933 

11.6167 

13.90 

10.192 

0.0 

1.0827 

2.1650 

3.2080 

0.3307 

5 . * 1 3* 

6.5 

7 • 6b  3 

16.0 

2 

0.0 

2. 3233 

0.6067 

6.9700 

9.2933 

11.90 

13.90 

1 * .2825 

0.0 

1.0827 

2. 1650 

3.2060 

0.330  7 

5,013* 

6.5 

8.079 

18. 

2 

0.0 

2.106 

0.29 

6.035 

8.58 

10.726 

12.87 

1 J.2005 

0.0 

.9996 

1.9991 

2.9987 

3.9563 

*.9979 

6.0 

7.095 

20.0 

2 

0.0 

1.7667 

3.5733 

5.3600 

7 . 10t>7 

8.933 

10.72 

1 0.9055 

0.0 

J 

.8326 

1.6652 

2.0976 

3.3300 

0. 1629 

5.0 

6.912 

c 

10O.0 

15.0 

0.2 

le.5 

1 1.0 

3.8 

.065 

76.0 

0.0 

0.2 

10.5 

1 1 .0 

1 1 .* 

• 065 

so.o 

15.0 

0.2 

9.75 

6.3 

2.  1 

.065 

76.0 

0.0 

0.0 

21  . J 

12.5 

12.5 

.12 

90.0 

6.0 

0.0 

13.0 

7.0 

6.0 

.12 

-9.75 

10.5 

2.9 

6 7.75 

5.0 

10. 0 

.5 

.05 

17.05 

1 .05 

17.*5 

0.0 

0.0 


The  first  three  pages  of  output  are  basically  a listing 
of  input  data.  On  the  next  page  are  the  principal  coef f icients 
of  the  roll  equation;  at  each  frequency  the  foil  coefficients 
form  the  first  line,  with  the  hull  coefficients  immediately 
below . 


C.  SAMPLE  OUTPUT 


A sample  case  of  FHROLL  output  is  given  below.  This 
output  results  from  the  above  input  data  and  is  fairly  self- 
explanatory.  Running  time  is  about  100  seconds  on  a CDC-6400. 


Sway,  roll  and  yaw  transfer  functions  are  then  listed, 
with  phases  relative  to  wave  elevation  at  the  CG.  The  final 
three  pages  give  root  mean  square  values  of  roll,  yaw,  flap 
angle  and  sway  in  the  three  specified  sea  states;  also  output 
are  absolute  motions  at  the  locations  specified.  The  quad- 
ratic regression  spectrum  of  Ref.  6,  obtained  by  analyzing  295 
wave  spectra  measured  at  station  'India'  in  the  North  Atlantic 
(59°N,  19°W),  is  used  in  the  irregular  sea  computations. 


N 

ii 


0 EL  HCG  XCG  HWG  YN6  OISW 

10.0000  ISO. 0000  2.0000  95.2500  .0860  .ifbnO  400.0000 

NFH  = IB  FW1  = .300  DFK=  .100 

N5E  A = 3 NWOS  = B 

hSW  TSh 

B.0000  8.6400 

10.0000  8.8700 

12.0000  9.1000 

XBOS  ZWOS 

( i)  0.0000  10.0000 

( 2)  S.0000  9.0000 

( 3)  10.0000  H . 0000 

( 4)  IS. 0000  7.0000 

( 5)  20.0000  6.0000 

( b)  10.0000  1.0000 

( 7)  10. 0000  15.0000 

( B)  10.0000  22.0000 

GM I N=  0.0000 


F I A V = -0.0000 

YAwAV=  -0 

.0000 

S«*  A Y A V=  -0. 

oooo 

STATION  2.00 

ABSCISSAS 
0.0000  .4592 

. 9 1 B4 

1.3776 

1.8369 

tt 

2.2961 

2.7553 

UHOINATES 
0.0000  .7379 

1 .4759 

2.2138 

2.9519 

3.3B9H 

4.427b 

STATION  4.00 

ABSCISSAS 
0.0000  .8639 

1.7279 

2 .59 1 B 

3.4558 

4.3197 

5.1 836 

OWD INATES 
0.0000  .75/4 

1.5949 

2.3922 

3.1897 

3.5871 

4.7H4S 

STATION  6.00 

ABSCISSAS 
0.0000  1.2770 

2.5540 

3.8310 

5.1080 

6.3850 

7 .662 0 

OBO INATES 
0.0000  .8569 

1.7137 

2.5706 

3.4275 

4.2843 

5.1412 

Who 
1 « 99011 


3.2146 
S • 1 6o0 

6.0476 

5.5b2<i 

8.9390 

S.99B0 


32 


STATION  «. 00 
AbSC I SS AS 


0.0000  1.6291 

3.2583 

4.8874 

6.5166 

8.1457 

9.7749 

1 1 .4040 

OHDINATtS 
0.0000  .9172 

1.8344 

2.7516 

3.6688 

4.5860 

5.5033 

6.4150 

STATION  10.00 

ABSCISSAS 
0.0000  2.2417 

4.4833 

6.7250 

8.96b7 

1 1 .208 J 

13.4500 

13.5015 

OHUINATES 
0.0000  1.1052 

2.2103 

3.3154 

4.4206 

5. 585 7 

6.6300 

6 . 8 J 1 0 

STATION  12.00 

ABSCISSAS 
0.0000  2.3233 

4.6467 

6.9700 

9.2933 

11.6167 

13.9400 

14.1020 

OKUlNATtS 
0.0000  1.0B27 

2.1654 

3.2480 

4.3307 

5.4134 

6.5000 

7.2470 

STATION  14.00 

AbSCISSAS 
0.0000  2.3233 

4.b467 

6.9700 

9.2533 

11.6167 

13.9400 

14.1920 

OHO  I NATES 
0.0000  1.0827 

2.1654 

3.2480 

4.3307 

5.41 J4 

6.5000 

7.6630 

STATION  lb. 00 

ABSCISSAS 
0.0000  2.3233 

4.6467 

6.97U0 

9.2933 

1 1 .9400 

13.9400 

14.2825 

OHOlNATtS 
0.0000  1 . 0U2  7 

2.1654 

3.2480 

4.3307 

5.41  JO 

6.5000 

6.0790 

STATION  IB. 00 

ABSCISSAS 
0.0000  ? . 1450 

4.2900 

6.4350 

8.5800 

10.7250 

12.8700 

1 J . 2 4 0 8 

UHOINATES 
0.0000  .9996 

1 .9991 

2.9987 

3.9983 

4.^979 

6.0000 

7.4950 

STATION  20.00 

ABSCISSAS 
0.0000  1.7867 

3.5733 

5.3600 

7.1467 

8.9330 

10.7200 

10.9456 

OHDIWATES 
0.0000  .8326 

1 .6652 

2.4978 

3.3304 

4.1629 

5.0000 

5.9120 

1 

1 
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XUCT  Mt-AN  bCoJAMfS  IN  SEA  S T A T t A 


SStA 


e.uo 


Tbt  A = A .6* 


, 


7llLL  A i^l ) YAW 
D I S2L  ACE  Mt  NT 
OtG 

VELOCITY 

UtO/btC 

ACCELERATION 

Ut6/SEC°*7 

MOLL 

2.379 

1 .6  74 

1.473 

Y Aar 

.39? 

.26$ 

.219 

flae 

3.322 

2.907 

3.824 

S«*AY  AT  POSITION 
D I SFL  ACE  MEN  T 
FT 

lNDlCATtl) 

VtLOCITY 

FT/SEC 

accelekat ion 
1 T/SEC*»2 

Cl) 

1 .**h 

.977 

.828 

Xs 

0.0 >1=1 0.0 

1 .404 

l .046 

.928 

Us 

3.0  • 7 = 9.  0 

1 .847 

1 .069 

.9  34 

X = 

10.0./=  h . o 

1 .628 

1.110 

• 9S9 

A: 

lb. 0.7  = 7.0 

1 . 733 

1 .169 

.948 

X = 

20.0.7=  6.o 

1 .«SH 

1 .24  3 

1.081 

X = 

10.0.7=  l.o 

l.*19 

.966 

.828 

X = 

1 1) . 0 . 7 = 1 8 . 0 

1 .8*39 

1 .271 

1.106 

X = 

10.0.7=72.0 

2 . 1 OS 

1 .44  3 

1 .262 
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hLuT  PF  AN  SUOAMtS  IN  SEA  S T « T t b 


HbE  A 


10.00 


TSEA 


H.a  7 


I 


POLL  AM)  YAn 
C 1 SPL  *Cr  ME  N T 
l)EG 

velocity 

uto/sec 

ACCtLEPA 1 ION 
OtG/SEC®*/ 

MOLL 

3.0S4 

/.  100 

1 .821 

1 

r An 

• SOM 

.33b 

.27  2 

1 i 
■ ; 

flap 

A.  187 

3.S88 

4.3S2 

I 

SnAY  ai  position 
ll  I SPl  ACEMF N T 
F T 

IMMCATtU 

VtLOCITY 

FT/SEC 

MCCELtPAT ION 
FT /SEC**/ 

LG 

1 • b H 3 

1 .237 

1 .02H 

A = 

o.o«/= io.o 

I .0  30 

1.317 

1.1*4 

j 

b . o « / * p.o 

/.00b 

1.3*0 

1.1S7 

J 

1 o . 0 » / = 8.0 

/.lib 

1 .*0b 

1 .IPO 

AS 

J 

• 

II 

■>/ 

• 

O 

2./S8 

1 .*H3 

1 .2*1 

A s 

j 

/0.0 */  = b.O 

/.*?* 

1 .S7V 

1.308 

A a 

10.0,/=  1.0 

1 .84* 

1 .222 

1.022 

i 

A = 

1 0 . 0 , / = I s . 0 

/.*lb 

1 . bOH 

1.373 

i ** 

10.0,  /■/’<*  a II 

/ . T 3/ 

1.82b 

1 .Shb 

HOOT  ME  AN  SQUARES  IN  SEA  STATE  5 


HSE  A 


12.00 


T St  A 


9.  10 


ROLL  AND  YAK 
displacement 
0E6 

VELOCITY 

OEb/SEC 

acCELEWaT ION 
DtO/SEC**? 

ROLL 

3.734 

2. SI  1 

2.140 

YAM 

.628 

.404 

.322 

FLAP 

4.984 

4.215 

S.072 

SMAY  AT  POSITION 
DISPLACEMENT 
FT 

INDICATED 

VELOCITY 

FT/SEC 

ACCELERATION 

FT/SEC««2 

CG 

2.338 

1 .497 

1.217 

x = 

0.0«Z=10.0 

2.378 

1.S82 

1.350 

X = 

5.0»Z=  9.0 

2.478 

1.624 

1.367 

x= 

10.0»Z=  8.0 

2.622 

1.696 

1 .407 

x = 

1 5 • 0 «Z  = 7.0 

2.802 

1.794 

1 .469 

X * 

20  • 0 »Z  = 6.0 

3.013 

1.913 

1 .550 

X* 

10.0.Z=  1.0 

2.286 

1.47S 

1 .209 

X: 

10.0.Z=15.0 

2.990 

1 .941 

1 .623 

X* 

10.0.Z=22.0 

3.379 

2.201 

1 .851 
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D. 


COMPUTER  PROGRAM  LISTING 


A 


| 

3 

i 

- 

j 


A complete  listing  for  FHROLL  follows.  It  is  worth 
noting  that  FHROLL  departs  slightly  from  Ref.  1 in  using  the 
methods  of  Jones*  to  calculate  C(k)  and  Se(k);  this  modifica- 
tion is  made  because  Jones'  formulation  takes  aspect  ratio  into 
account.  Another  noteworthy  point  is  that  in  calculating 
strut  roll  damping  terms,  account  is  taken  of  the  variation 
in  roll  velocity  along  the  strut's  span. 


* Jones , R.  T.:  The  Unsteady  Lift  of  a Wing  of  Finite  Aspect 

Ratio.  NACA  Report  681,  1940. 
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PROGRAM  FhHOLL  ( INPUT, OUTPUT. T AH t 5= INPUT, TAHE2= INPUT. TAPES =OUTPUT> 
COMMON/COM  1 /OO.G ( 10. 1 1 ) 

COMMON/NEW/XA (25) .OXA12S) . XCG  ,tL . NST , MCG ,CA Ah ,0  ISP , MmO 
C0MW0N/NEW2/W.U.A22.B22.A2A.H2A , a26 .826  . AAA  ,H  A A . A Ab  .Ha.,  . 
lA62,862»A6A«bbA,A66»866»EF  (10) 

C0MM0N/NEW3/A22F.A2AF,A?6F  .AAAF.AA6F.A66F 

COMMON/NEWS/OAM (6)  ,S(6)  ,SX(6)  »T  (6)  «Z(6)  , NFS  » d ( 6 ) .COSS  <6)  .SINU  (6) 
C0MM0N/NEW6/FIAV,  YAWAV.SwAYAV 

COMPLEX  A I .d22.b2A,b26.C2A,C26.F2.UAA,sA6.CAA,CAh,FA.t)62.tJ6A,rtb6. 
lC6A.Cbb,F6,CK.SE.UI .30.H2F ,C2F .UaF .CAF.U6F , C6F .ri2H « C2H . RAH  . 

2CAR .HbN.C6H.Alw 

COMMON  PI .HP  I .UP  I ,TPI .MD.MODE.OPh.CR.RAT.SUR.OfcG. I S T , DHT ,HPh. SG . N 
10E, POM, VOL. Uc*i, UN. OMEGA, CP.wVh.lt),  OOG  «IG,XX(25.7)  . Y Y ( 25 . 7 ) .DEL  ( 25  . 
27) . SNE (25.7) ,CSE<25.7) ,FP(7) ,dL 00 (25,7,7) . YLOG (26,7.7) .C0M1A.1)  «c 
3T(1A.1a).PSI1(7.7).PSI2(7.7).PHA(7).PHV(7) 

DIMENSION  XP0S( 10) «ZPOS (10). SwEEP (6) , ALF (6)  ,CH (6 > . TC  < 6 ) . 

ICE (6) , Cd  AH ( b ) • X ( 6 ) .H(b) « CL A ( b ) .CLH(6) .A0(6),ZPl(b).HB(6>, 

2A(6) » SC ( 6 I » YZ (6 ) .OUTM(AO.IO) .OSP( 15) , 

3VEL (15)  .ACC (15) .SPEC! 10) ,SPP (2) ,YA(S) ,HSW l 10) ,TSw( 10) 

PI=3.1A1 5927 
TPI =2 . *P I 

>9  HEAD  13.U,EL»HC6,XCG.RHG.YhG,DISP,RHO 
IF (EOF ( 5L INPUT) .NE. 0.0) STOP  1111 
WHITE  101 

WHITE  13.U.EL. HCu , XCG . PHG .YhG.DISP .WHO 

READ  A0.NFH.FH1 ,DFH 

WHITE  A3.NFH.FR1 ,DFR 

READ  50.NSEA.NPOS 

WHITE  51.NSEA.NH0S 

IF(NSEA.LE.O)GO  TO  67 

WRITE  53 

00  52  1=1 .NSEA 

READ  13.HSWII) .TSw(I) 

52  WRITE  13.HSWII) , T SW ( I ) 

IF(NPOS.LE.O)  GO  TO  67 
WRITE  1009 
DO  82  1=1. NPOS 
READ  13.XP0SII) »ZPOS(I) 

12  WHITE  lOOl.I.XPoS(I) ,ZPOS(I) 

37  CONTINUE 
CALL  HULL  I 
NFOIL  = 3 
READ  18.NSTRUT 
WRITE  101A, NFOIL. NSTRUT 
NFS=NFOIL*NSTHUT 
WRITE  1002 
00  83  1=1, NFS 

HEAD  13. GAM ( I ) .SWEEP ( I ) .ALF ( I ) .d ( I ) .CR ( I ) ,CE ( I ) »TC ( I ) 

WHITE  13.GAMI I ) .SwEEP ( I ) . ALF ( I ) ,a ( I ) ,CR ( 1 1 ,CE ( I ) . TC ( I ) 

HEAD  13,X(A) ,Y(a) ,m(a> ,X(5) .H(S) 

WHITE  1003. X(A) .Y(A) , H ( A ) *X(5) »H(5) 

IF  (NSTRUT .LE • 2) OU  TO  2A 
X (NFS) =X (A) 

Y (NFS)  = 0.0 
H (NFS  ) =“.*j*tJ  ( NFS ) 

?A  CONTINUE 

READ  13.  ttF.PF.EFF.wF.ZETF 


fl  ' 1 
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Wk  I TF  1 0 0 4 • h(-  « MF  « H.F  F • l*»F~  • /F  1 F 

H£  AO  1 .1  , NM  . /(_  T rt 

rthITc.  100b»-»w«/F.  Th 

kfc'AO  l j.UF'UD*GFU,(jF 

Wf<  I TF.  I 0 06 . UFOU • OF  0 « OF 

Kt  AO  ljt.QkUO.'iHD.  Uk 

wklTK  IOO/.OnDO.UmO.OK 

YF  = n l 1 ) -r*K 

CF  = C« ( 1) - (CM ( 1 ) -CF. (1) ) » ( YF ♦.S«cK) /rf ( 1 ) 
koT  OiHKOkML  ANGLES  IN  CGNVHNT  10  <AL  FOHM 
G1  = («AM(1)-Gu.)/S7.3 
GAM  ( 1 > = IHd.-oA M(  I ) -Gam  (4) 

GAM ( ?) = GAM (2) -6AM (4) 

GAM  CD  =O0.-GAM( j) 

G A *M  ) = - 1 ? A M ( 4 ) 

GAM  (S  ) = -H0. 

GAM  ( 6 ) =-S»0  . 

ChANC-fc  ANC-LtS  FkOM  OEGwptS  TO  haiJIANS 

00  l 1=1. NFS 

SwF.fck  ( I ) -SWtt.M  I I)  /67..T 

k = G A M ( l ) / 7 . J 

SING ( I ) =S IN  (k ) 

COSS(  I) =CGS (k) 

DO  2 1 = l « NF  s 

ALF  ( I ) =AI_F  ( 1 ) /s7 . 2 

CHAK ( I ) =.S"  (CH  ( I)  tCt  ( I ) ) 

COM  I MJt 

DU  T 1 = 1.  NFS 

S ( l ) =0  ( I ) »CDAK  ( I ) 

H INT*H ( A) *ri ( A) *«HS (SING (») ) 

F(<fW  = k<6)  *hI*IT*b<P)  *SlNG<2> 

7F  =-kCG-h  I NT  ♦ ( YF  ♦ .4<>»HF)  »SlN0  ( 1 ) 

YF=Y  (4)  ♦ ( YF*.4  2°HM  ”C0SS  ( l ) 

0=  (H  ( l ) y>COSS  ( 1 ) *r>.  (?)  SCGSS  (2)  ) *«<•/ (S  ( 1 ) “GOSS  ( 1 ) *S  ( ? ) *COSS  ( 2 ) ) 
T = b ( 2 ) oCOSS  ( e7 ) — Y (4) 

IF  (AHS(T)  .Lt.O.ODSO  TO  4 
I SkL I T = 1 

H ( = h INT  ♦ *SIN«  (2) 

' <2>  = Y (4) -,42"rt (2) *COSS (2) 

A ( 2 > * U 
GO  TO  s 

iski.i  r = o 

H ( 2 ) = Y ( '*  ) /C0SSI2) 

H (2)  =H  I NT  ».S*rl  (2)  ASIi.G  (?) 

Y (?)  =Y  ( '.)  - ,S*H  ( 2)  *COSS  (2; 

A (2)  s2.*<) 

CONTINUE 

Fi  ( I)  =nIM-.  42*6  < 1 ) USING  ( 1 ) 

H(  J)  at,  ( •>)  -.  A2»tl  (!)  <>l,ING  ( .1)  *H(S) 

H ( A ) =1I\T-.S»M («*) «AkS(SlM,  ( 6) ) 

k ( S I = . D 0 » 0 ( *.> ) » H ( ‘j ) 

X ( | > = A ( ? ) = X l A | 

X ( J ) a X ( ) ) 

00  m 1=1. NFS 

SX  ( I ) -X  ( I)  ( I ) 

/(!)=-.  it(--kU) 

(ll|  ■:  Y (A)  » .<*?  <‘M  ( l ) CUSS  ( l ) 
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Y(3)=.42*8<3)*C0SS<3> 

Y (4)  =Y  (4)  -.50*8  (4)  »COSS  (4) 

Y<5)=0.0 
A ( 1 ) = A<2) 

00  7 1=3, NFS 
7 A(I)=2«*b(I)/CbAH(I) 

EM=2240.»OISP/ <32.2*HH0) 

U=1 .689*U 
XI  = (RRU*EL> **2*£M 
21=  ( YRG*EL>  **2*EM 
AI=(0. 0,1.0) 

C LIFT  CURVE  SLOPE  CALCULATIONS 

C 

C INCLUDE  FREE  SURFACE  EFFECTS  FCR  FOILS  I ANO  3 

00  22  I = 1 »NFS 
CL A ( I ) = CLH<  I)  = 0.0 

22  AO  ( I ) = TPI*< l.-.96*TC(I> ) *COS ( SWEEP ( I ) ) 

IF  (U  .LE.  0.0)  GO  TO  28 

00  23  I = 1.3,2 

HCSQ  = 20.0*(H(I)/C8a*II) 1**2 

AO  ( I ) = A0(H*(I.*HCSQ)/(2.*HCSO) 

HS  = H<I)/3(I) 

2P1 ( I ) =1 .0*8IPL (HS) 

23  CONTINUE 

2P1  (2)  = ZP1 (5)  = 1.0 
C STRUT  END  PLATE  EFFECTS 

HB(4)  = 1.9*b<l)*C0S<Gl)/B(4) 

HB ( 5 ) = 1.9*8(3)®C0SS(3)/B(5) 

IF  (NSTRUT  .LE.  2)  GO  TO  25 

HB<6)  = 1 .9*  ( till ) *COSS  ( 1 I *b  <2)«C0SS  (2))  /B  (6) 

ZP 1(6)  = 1.0 

25  CONTINUE 

00  26  I = 4, NFS 

26  A ( I ) = A ( I ) * ( 1 .0*M8 ( I ) ) 

HS  = Y (4) /B (4) *2** (3-NSTRUT) 

ZP 1(4)  = 1.0*bIPL(HS) 

00  27  I = 1 .NFS 

CL A ( I ) = CLALFI AO ( I) , A(  I) ,ZP1  ( I) ) 

27  CONTINUE 

00  9 I = 1,3.2 

HCSQ  = 20.0*(1.05*H(I)/CBAR(I) ) **2 

AO ( I ) = TPI*«1.-.96*TC<I) )*C0S(S«E£P(I))*(1.*HCSQ)/(2.*HCS0) 
HS  = 1 • 05*H (I)/B(I) 

ZP1 (I)=1.0*3IPL(HS) 

CLH(I)  = ALF  ( I ) * (CLALF  (AO ( I) . A( I) ,ZP1 ( I) )-CLA(I) ) /{ .05*H( I) ) 
9 CONTINUE 

20  CONTINUE 

00,91  1=1, NFS 
CLH(I)=S(I)*CLH(I) 

SC  ( I ) =5 ( I ) *CL  A ( I ) 

YZ(I)sYM)  »COSS  ( I ) *1  < I ) *S  I Mi  ( I ) 

91  CONTINUE 

00  92  1=4, NFS 

YR  = Y ( I ) -.5*8 ( l ) *CCSS ( I ) 

ZR=Z< I) *.5*8 ( I) »ABS(SING( I) ) 

YHH  = YH*COSS ( I I ♦ ZH“S I NG ( I ) 

Y4(I)=((YRR*B(I)) ** J-YRH*»3 ) / ( 3 . 0*H ( I ) ) 
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92  CONTINUE 

DO  93  1=1.3 
43  Y 4 ( 1 ) =Y t ( I ) **2 
d (5) = .S*0 (5) 

SC(S)=.B»SC(5) 

IF (NSTRUT.LE.21GO  TO  130 
H (NFS) = .5»B (NFb) 

SC (NFS) =.5*SC(NFS> 

130  CONTINUE 

C COMPUTt  FMtOUENCY  INDEPENDENT  TERMS 

A22=A24=A26=A44=A46=a66=0.0 
00  10  1=1, NFS 

0=TPI*b( I ) * (CdAR< I ) /2.) **2 
00  = 0* S I NO ( I ) **2 
A22=A22.0Q 

A24  = A24-Q*S ING  < I) »YZ( I ) 

A26=A26.00*SX(I) 

A44=A44*0*Y4 ( I ) 

A46=A46-Q*SX ( I ) *S I NG ( I ) * YZ  ( I ) 

A66=A66*00«SX ( I ) **2  »TP I *B ( I ) *CB  AM ( I ) **4/ 1 28 . *S I NO ( I ) «*2 
10  CONTINUE 
A64=A46 
A62=A2b 
A22F=A22 
A24F=A24 
A26F  = A26 
A44F  = A44 
A46F  = A4G 

A66F  = A66 

YZF  =YF *COSS  ( 1 ) *ZF*SING( 1) 

CALL  FLAP/RFfTl,T4.T7.T0,TlO,Tll) 

A2F=-2.*BF*T1* (CF/2. ) «*3 
A4F  = - A2F* Y ZF 
A2F=A2F*SING(1 ) 

A6F  = A2F*SX ( 1 ) *2  • *0F  * (CF/2. ) **4» ( T 7 *PF*T 1 ) *SlNG( 1 ) 
A2R=-TPI*8 (5) *C0AH (5) **3/ lb. 

A4R=-A2H*Z (5) 

A6H=A2R*SX (S) *TPI*B (5) *CBAH (5) **4/ 128. 

WRITE  1042 
WHITE  1043 
00  99  IFH=1,NFK 
W=FR1 ♦ ( IFH-1 ) *OFR 
OW=W*w/32 • 2 

B22=824=C24=B2b=C?b=F2=<0.0,0.0) 

S44=C44=U46=C4b=F4= (O.O.O.O) 

062=0b4=C64=0b6=Cbfi=F6= (O.O.O.O) 

88=B (4) 

00  20  1=4,5 

IF ( I . tO . 4 ) GO  TO  19 

8B=2.*ri  (5) 

19  CONTINUE 

Q = P 1*8(1) »C8AM( I ) «SDAPP (w.HH) 

822=B22*0 
844=B44*0*Y4 ( I ) 

H6b=Bbb*Q*SX  ( I ) **2 

20  CONTINUE 

IF  (NSTMUT.Lt.2IG0  TO  133 


44 


1=NFS 

'J=*»I*FHl)«CdAN(l):*'»UA»«F'<*»~ttt» 

H22  = H<;2  »U 
044=R44*(JOY4  l I ) 

B6b=B5o*OoSA  ( I ) **? 

133  CON T [NluC 

IF  tU.liT  .0.0)  00  TO  16 
CALL  ZE  NO ( rt * 322 . b4 4 • b66 ) 

FJ2F  = R4F=H6F  = C2F  = C4F=C6F=(0.0.U.U) 
b2R  = B4K  =36Ft  = C 2F<  = C4H  =C69=  ( 0. 0 . U . 0 > 
GO  TO  17 
16  CONTI NUF 

DO  11  1=1. NFS. 

0=.5aCBAR ( I ) aA/U 

CALL  THE0J0N<A< I ) .G.CK.SE) 

F>  = SING  ( I ) '•>»,? 

H = SING< I ) *Y Z ( I ) 

T = SX ( I ) -CHAN ( 1 ) /a. 

Q I=0*SC ( I ) ttCK 
U0=F*(1I 
B22=B22*OU 
B62  = Hb2*OU‘»X  ( t ) 

B24=324-OI»N 
b44=B44+GI«Y4 ( I ) 

H64=Bb4-<) r-*H»X  l I ) 

HA  = -0»Tf,I  «H  (1)  '*  (CHAN1  I)  /?  .)  A*  2 
0 U = 0 I a T 

N26  = H26  ♦ F*'*  (dA*iJIJ> 

B46=B46-Ra (BA.OLl 
B66=B6fa*B° ( Ta6A*X ( I ) ° 0 0 ) 
QI=-U°a2aCLH( I)«CK»Y(  1) 

C2A  = C24.0IaSlNu(  1) 

C44=C44-QIaYZ  ( 1 ) 

C64=C64*llI»SING  < 1 ) »X  ( I ) 
OI=-Uoa2aSC ( I ) »CK«SlNG ( I ) 
C26=C26*0I»SING ( I ) 

C46=C46-QIaYZ  ( I ) 

C66  = C66*0IaSINb ( I ) *X  ( I ) 

IF  (I  . LE  . 3)  GO  TO  94 
IF  (I  .GT , 4)  GO  TO  95 
YN  = Y(I)-.5aB(I)aC0SS(I) 

ZH  = Z ( I ) * .5*3 ( I ) *AHS (SING(  I ) ) 

YHK  = YR*C0bS(l) ./N»SING( I) 

BH  = h l l > 

HO  = h(I)-.5*B(I>°'>BS(  SING  ( I ) ) 

HT  = h( I) *.5*d< I)«AbS(SING(  I)  ) 

GO  TO  96 

95  CONTINUE 

FtH  = 2. 0.1(1) 

YRH  = -Z ( 1 ) - 3(1) 

HO  = H ( I ) - Bit) 

HT  = HO  ♦ By 

96  CONTINUE 

IJWW  = 0»o  AbS  ( G I I.G  I I ) ) 

YTT  = YNN  ♦ 3y 
Tft  = (jww*YRn 
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tt  = uwwoytt 

T = EXP  ( TR-Q«*H() ) / (ort«UWW**?) 

P = -T» (EXP (-TR) *( TR* 1 . ) -EXP (-TT) * (TT.  1 . J ) 

R = -SING(I) 

T = -Y  < I ) 

TT  = EXP(-OW«HO)-EXP(-Gw®HT) 

TT  = TT/(Qw»<HT-hO)  ) 

DO  97  J=  1 ♦ 2 
P a -P 
R = -R 
T = -T 

01  = -.b«U«CtXP  ( AI»Ow«T>  *SC  ( I ) (R*AI*COSS  ( I ) 1 

F2  = F2  ♦ QI»R»TT 
F4  = F<*  - GI*P 
F6  = Ft  * QI*H*X < I) *TT 
97  CONTINUE 

GO  TO  98 
9*  CONTINUE 

P = -YZ < I ) 

R = -S  I NG ( I ) 

T = -Y  ( I ) 

DO  12  J=1.2 
P = -P 
R = -R 
T = -T 

QIs-.S^U^CExP  («  [ ® G w » T ) * ( U ® C L H ( I ) »CK»SC(  I ) »Sfc**i<»EXP  ( -Qw«H  ( I ) ) * (R 
1*  M»COSS  ( I ) ) ) 

F2=F2*0I*R 
F 4=F4-G  I *P 
F6  = F6*0I»H»X ( I ) 

12  CONTINUE 
98  CONTINUE 

IF ( I .NE. 1 >G0  TO  11 
P = YF®COSS ( 1 1 ♦ZF®SINb( 1 ) 

Q I = T 4 -CL  A (1>»CK*T11/TPI 
ui=-.s®u»eF*cF®»2®oi 
02F=OI»SING( 1 ) 

B4F=-GI *P 

GI=T4*bX(  ll-CLAl  1 ) “ CK  6 T l 1 ® X ( ll/TPI 
0 I =01  ♦ ,5“CF  ® ( !i-T«.PF»H..5Mlll 
BftF=-.5«U«SF*CF®»2®UI®SING( 1 ) 

QI=U**2*0F»CF»CLA( 1)*CK*EFF 
C2F=0I»SING'1> 

C4F=-GI *P 

C6F=C2F*X  ( l ) - .b»U«»2*BF»CF**2«  (1  *.♦  T 101  »SING  ( 1 ) 

11  CONTINUE 

0 I = C L A ( 5 ) *CK 
P = -.25*U*S<b>  *C»AKIb) 

B2R=P® (“I *01 ) 

B4R=-i2P»/ ( S ) 

B6W=P‘>  (PI®  (bX  (b)  - .25®  CRAW  (SI)  *.I*x  I SI) 

C?K=-.b®U»«2®b (S) *01 
CARa-C2P»/«b) 

C6H=C2R»x (b) 

17  CONTINUE 

WRITE  l!)40i</tA»i»l  <H<i‘t»C*|t*l  * 

C44=C*4*C»4H 


46 


CALL  HULL  /i 

L CALCUL  aTICiV  Or  A • 5 . ri  • b . L • .4  4iyu  r • n NO  I CUMML.  ir 

L COMMOTE  rtYGHOUYNAM  1C  MATRIX 

A I « = A I » * 

2 = * » w 

GU=-«2« < A22*EM)  »AI**rt22 
CALL  MATG(l.l) 

GU  = -w2:>«24*4Iw<>rt24*C24 
CALL  MaTG(I.J) 

Q0  = -w2<*A26*AI**U2b*C26 
Call  matgii.si 
OQ  = -w2«A2K*AIW»rf2F«-C2F 
CALL  MATGU.7) 

OU=-W2*A2H»A  I*»>$2h*C2R 
CALL  MATG (1*9) 

GQ  = -W2aA24*A[W»rl24 
CALL  MATGI3.1) 

GU  = -*2<MA44*XI>  *A Iw*SA4*C4m 
CALL  MATG(3*3> 

QG=-42uA4b*A  I H»h‘lh*C4h 
CALL  MATG(3,b) 

Q0  = -w2*A4F*AI»*tl4F*C4F 
CALL  MA  TG ( 3 * 7 ) 
gG  = -W20A4H*AIW*ti4H*C44 
CALL  MATGI3.4) 

QQ  = -4  2°Ab2*AI<l“Bb2 
CALL  M A T G ( 5 * 1 ) 

0a  = -/»2»Ab4»AIw*ribA*Cb4 
CALL  M A TG ( 3 • 3 ) 

O0  = -V«?°  ( A6lW  I ) *A  IW*b66*Cbb 
CALL  w A TG ( 5 * S ) 

UU=-'*2*AbF*  A I **HbF  *CoF 
CALL  M A T G ( 5 * 7 ) 

0U=-W2*A6M»A I *»RbR*C6H 
CALL  MATG(3*4) 

OU=  (-W2*uFUD*  AlW*OFO*OF  ) « <-*nF*«2) 

CALL  MATGI7.3) 

OQ=-w2«AIW«2.*ZtTF*AF*wF*A2 
CALL  M A TG  ( 7 * 7 ) 

OQ  = ( -W2*OHOO*  A 1 w«OHD*OW) * (-«K**2) 

CALL  MAT G ( 9 « S ) 

OQ=-W2*  A I w<*2.  */ETR*WH*wk**2 
CALL  MATG(9*9) 

OG= (0.0 .0.0) 

CALL  M A TG ( 7 * 1 ) 

CALL  M A TG ( 7 *5 ) 

CALL  Mm  TG ( 7 * 0 ) 

CALL  MATG  (4*1) 

CALL  M A TG ( 4 ♦ 3 ) 

CALI.  MA  TG  ( 4 * 7 ) 

C COMMUTE  EXCITING  FORCE  VECTOR 

EF  < l ) =rtE  AL  ( K2 ) * ji F ( 1 ) 

EF (2) =AIMAG (F2) *EK (A) 

EF (3) = RE  AL (FA) *EF (3) 

EF ( A) =AlMAG (Fa) »eF (A) 

EF  < b ) =KEAL  (Kb) *eF  (b) 

EF  (6)  =A1MAG  (Kb)  *EF  < r> ) 
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I 


E 


DO  14  1=7.10 

14  EF ( I ) =0 .0 

C SOLVE  FOR  MOTIONS 

DO  15  1=1.10 

15  0(1.11) =EF ( I ) 

WHITE  104l.A44.B44.L44.tF (3) .EF (4) 

CALL  SOLV (G.EF.10. 1 l.INDX.ICK) 

DO  225  J=1.10 
225  OUTM ( IFH. J) =EF ( J) 

99  CONTINUE 

C OUTPUT  FHEQUENCY  RESPONSE 

WRITE (6.212) 

WRITE(6«231) 

WHITE (6.214) 

WRITE (6.215) 

DO  227  L W = 1 ,NFH 
W=FR 1 ♦ (Lw-1 ) *DFR 
WL=TPI*32.2/W**2 
WSLP=TPI/WL 
WL=WL/EL 

SAMP=SQRT (OUTM(LW, 1 ) »*2 *OUT M ( L w . 2 ) »*? ) 
SPH=57.3«ATAN2(OUTM (LW.2)  .OuTm(Lw.I) ) 
HAMP=SQHT (0UTM(LW,3)**2*0uTM(LW.4)»<»2>»B7.3 
RPH=57.3*ATAN2(0UTM(LW,4) ,0UTM(LW,3) ) 

YAMP  = SUHT  (OUTM  (LW  ,5)  i*»2*0UTM  (l*  .6)  <**2)  »57 .3 
YPH=57.3“ATAN2(0UTM(Lw,6).0UTM(Lw,5)) 

WHITE (6.216) W.SAMP.SPH.H AMP,  RPH.Y AMP, YPH.WL 
227  CONTINUE 

IF (NSEA.Lt.O) GO  TO  1000 
DO  54  JS= 1 »NSE A 
H5EA=H5W ( JS> 

TSEA=TSW(JS) 

IF (HSEA.GT.0.0)  GO  TO  30 

ISEA=0 

GO  TO  35 

30  IF (HSEA.GT.1.0)  GO  TO  31 
ISE A=  1 

GO  TO  35 

31  IF (HSEA.GT.3.0)  GO  TO  32 
ISEA=2 

GO  TO  35 

32  IF (HSEA.GT .5.0)  GO  TO  33 
ISEA=3 

GO  TO  35 

33  IF (HSEA.GT. 8.0)  GO  TO  34 
ISEA=4 

GO  TO  35 

34  IF (HSEA.GT. 12.0)  GO  TO  46 
I SE  A = 5 

GO  TO  35 

4h  IF (HSEA.GT. 20.0)  GO  TO  47 
I SE A=6 
GO  TO  35 

47  IF (HSEA.GT. 40.0)  GO  TO  4M 
ISEA=7 

GO  TO  35 

48  I SE  A = B 


! 
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I 


I 


I « 
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Is  CONTI  NUt 

ntohnpos.  i 

00  70  1 = 1.  MOT 
7"  USP ( I ) = VEL ( I ) =ACC 1 1 ) =0 
Ht  T D = bc.  T V = HE  TA=U.O 
00  71  LWsl.NFH 
w=FRl ♦ (LW-1  ) aUFW 

w4=w2»i*2 

FW=SE AST  (hSEA  . TSE  A . m) 

0=Fw» (OUTM  (Lai.  7 ) »*?*OUTM  (L-  ,8)  *»2) 

bETD=bETD.O 

betv=betv*u‘>w2 

BET  a = HET  A.G»w4 
00  72  1=1.3 
J = 2*I 

72  SPEC  <I)=F**(0uTm<LW,J-1)°*2»0iJTM(LW,J)»*?) 

00  73  1=1.3 

DSP ( I ) =OSP ( 1 I .SPEC  ( I ) 

VEL ( I ) = VEL ( I ) .SPEC  ( I ) **2 

73  ACC  ( I ) =ACC  ( I > .SPEC  ( I ) »W<. 

IF (NPOS.LE.O)  00  TO  71 
00  76  1=1. NPOS 

DO  77  J=  1 . 2 

7 7 SPP  ( J)  = 0uT8  ILw  , J)  -ZPOS  ( I ) »>R)T.m  (LW.  J»?)  - ( APCS(  I ) =EL/20  . CG ) »OUTP<L 

«W»J»4) 

7 6 SPF.C  ( I ) =Fw»  (SPP  ( 1 ) *°2»SPP (2) ««2> 

00  7 H J=4.NT0T 
OSP(J)=OSP(J) » SPEC ( J-3 ) 

VEL ( J ) = VEL ( J ) »sPtC ( J-3) *W2 
78  ACC (J) =ACC (J) .SPEC (J-3) «*4 

71  CONTINUE 

00  7*  1=2,3 

DSP  ( I ) =SQ*T  (0Fw»'J5P  ( I ) ) *S7.3 
VEL I f ’ = SQHT (OFPOVEL ( I ) ) *57.3 

74  ACCiO=.iQrfT(0FW*ACC(l)  )*57.3 
wHITE(6,2l7) ISEA.MSEA.TSEA 
WHITE (6.218) 

WHITE (6.219) 

WHITE (8.220) 

WHITE (6.221 ) USP (?) .VEL (2) .ACC (2) 

WHITE (6.22?) OSP(J) » VEL ( 3) .ACC (3) 

BETO=SURT (UFH«bETO) »B7.3 
HETV  = SOHT (OFH»bETv>  «57.3 
BETA  = SOHT (OFH»at  TA) »S7.3 
WHITE (6, 226  HE  TO. HE  TV. BET A 
OSP ( 3) =USP ( 1 ) 

VEL (3) =VEL ( 1 ) 

ACC (3) =ACC ( 1 ) 

00  7*>  I =3,NT0T 

OSP(  I ) =SOHT  (L)FH°OSP  ( I ) ) 

VEL ( I ) = SUH  T ( OF  8 « V f L (I)) 

7 S ACC ( I ) =SUHT (UFPOACC ( 1 ) ) 

WHITE (6.2?3) 

WHITE (6.219) 

WHITE (6.224) 

WHI Tt (6,226) OSP (J) .VEL (31 >ACC (J) 
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IF (NPOS.LE.OIGO  TO  54 

WRITE (6*240 M XROS  < 1-3) .ZPOS ( l-J) .DSP< I ) * VEL  ( 1 ) .mLC  ( [ ) . 1 =<♦  , NT  ,i  n 
54  CONTINUE 

1000  IF  (EOF  (5LINPUT)  ) SOSES') 

909  STOP 

18  FORMA  T ( I 1 ) 

212  FORMA! ( 1 H 1 / /25X *F HE QUtNC Y HESPONSE*) 

214  FORMAT  (//  15X4HSwAY  . 1 9X  4HR0LL  . c U A JH  Y A W .15XhR<*.L./L) 

215  FORMAT  (3X.  lrU,7X,3RAMP,5X,5HPHA6t. * 1 0 X . .in  AMP  . 5X  , 5hPm  aSE  * 1 0 < . 3"  AMR  . 
IX* 5HPh ASE ) 

216  FORMAT <F7.3*2F9.3*5X.2F9.3.5X,2F9.3.F15.3) 

217  FORMAT  ( im//10A*HOOT  MEAN  SQUARES  IN  SEA  S 1 A T E * 1 2 . 5X*hSE  A =»rS.2.5 
1X*TSEA  =*F5.2) 

2 1 B FORMAT  (///15X»HULL  AND  YA'**| 

219  FORMAT  t 15X»0ISPLACEMENT"1  lX“VtUOC  I TY'>  1 lX*ACCELt  RAT  IONJ) 

220  FORMAT (20X»0EG"1 5X*DEG/StC" 12X1 0RUEG/SEC»*2) 

221  FORMAT  1/6  X*HOLLuF  15.3  *2F2  0.3) 

222  FORMAT (/6X«YAw»F 16.3*2F20.3> 

223  FORMAT (///15X*SwAY  AT  POSITION  INDICATED*) 

224  FORMAT (20X"FT*16X»FT/SEC»13X9hFT/5EC»<*2) 

226  FORMAT  (/6X»FLAP*F15.3*2F20.3) 

228  FORMAT ( /5X<*CG*3XF  l 5 . 3 . 2F20 . 3 > 

231  FORMAT <//SX*S*AY  AMP  IS  NON-D I M E N S I ON AL * ROLL  AND  Y A * AMOS  IN  DEG/ 
FF T* ) 

240  FORMAT (/2X,*X=*F5. 1*.Z=*F4. 1 .F9.3.2F20.3) 

13  FORMAT (8F10. 4) 

40  FORMAT ( I2.2F10.3) 

43  FORMAT (/3X*NFR=*I3.5X«FR1=*F5. J.5X*DFR=*F5.3) 

81  FORMAT (2F10.4. 12) 

10  l FORMAT  < IHl  ,5X»U<>9X<*EL«8X*HCG*7X<‘XCG*6X0RRG«  7X* YRU*  7X *0  I SP»  7 X *RnO*  I 

1001  FORMAT </3X.» (*I2») °2F 10.4 ) 

1002  FORMAT  </5X*GAM*7X«SwEEP*5A*ALF»7X<»R»9X*CH»6A<HE*8X*TC») 

1003  FORMAT  ( / 1 X°X  ( 4 ) = *F8.4,5X»Y (4 ) = *F 7 . 4 *5X*m ( 4 ) = *F 7 . 4 . SX »X [ 5 ) = 
4 *F8.4*5X*H (5)  = *F7.4) 

1004  FORMAT  ( / 1 X *fc)F  = «F8.4*5X«PF  = «Fb.4*5X»EF  = »F6.4,5X»wF  = *FB.4, 
45X*2ETF  = »F7.4) 

1005  FORMAT  ( / 1 X*wR  = *F8.4.5X*ZETR  = *F7.4) 

1006  FORMAT  </lX»QFOD  = *F8 .4  *5X*0FD  = <*F8.4*5X*UF  = *F8.4) 

1007  FORMAT  </lX»UHDD  = *F8.4.5X*QRU  = *F8.4.5X*QR  = *F8.4) 

1008  FORMAT (/3X*HStA=*F 12.2.3X»TSEA=*F 12.?.3X,»NP0S=*. 13) 

10  09  FORMAT  (/13X,*XP(JS»10X*ZPOS*) 

1010  FORMAT  </3X*GAM ( I ) »2X*SWEEP ( I ) » . 4X . ° ALF ( I ) »6X . «M ( I ) *5X»  CR  ( t ) * 
15X,*CT(I) )*) 

1014  FORMAT  (/1X**NF0IL=  * I l * 4X  *NSTRuT  = * , II ) 

1040  FOHMAT  (F10.3.7E12.4) 

1041  FOHMAT  ( 1 0 X * 7E 12.4) 

1042  FORMAT  ( 1H1//,6X*R()LL  COEFF I C I ENTS » ) 

1043  FORMAT (//6X"W“10X*A44*HX«844R»HX»544l"8X«C44R»clX"C44l*9X>F4R»9X 
1 *F4 I */ ) 

50  FORMAT (213) 

51  FORMAT l/5X»NSEA=»I2,5X*NPOS=*I2) 

53  FOHMAT  ( /4X*H5V.*7X*>TSW*) 

END 
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SUdROU  T i NE  SULV (A.X.N.M, INOX. 1CK) 

c 

C SOLUTIONS  OF  N L INEAK  EQUATIONS  IN  N UNKNOWNS. 

C 

C H=N*i 

C 

C MATRIX  tQUATION  SOLVED  IS 

C 

C 8*Y=C 

C 

C WHERE  6 ( I . J=A ( I , J)  I *Js 1 *N 

C Y ( I ) =X ( I ) 1=1. N 

C C(I)=AU.M)  1 = 1. N 

C 

C IF  NO  SOLUTION  FOUND  IOK  IS  SET  EQUAL  TO  1 FOR  RETURN. 
C 

DIMENSION  A(N.M) »X(N)  , INOX ( N ) 

ICK  = Q 

OHIO  1=1. N 
INDX ( I ) =0 
10  X(I)=U.O 

00  20  J= 1 * N 
ZZ=  1 . OE- 1 0 
IROwsO 
DO  30  1=1. N 

IF  (INDX ( I ) .NE.O)  GO  TO  30 
TEST=AHS(A<I,J) » 

IF  (TEST .LE.ZZ)  GO  TO  30 

ZZ=TEST 

IROw=l 

30  CONTINUE 

IF  ( IRO« .tU  .0 ) GO  TO  20 
«0  INDX(IR0m)=J 
ZN= A ( IROM , J) 

II  =N*  1 

DO  50  K= 1 . 1 1 

SO  A < IROW.K) =A( IxGw.K) /ZN 
00  60  1=1. N 
IF  ( I .EQ . I ROW ) GO  TO  60 
Il=J*l 
1 1 =N* 1 

00  61  K = I 1 * 1 1 

A(I.K)=A(I.K)-A(I,J)»A(IROW.K) 

61  CONTINUE 
60  CONTINUE 
20  CONTINUE 

00  80  1 = 1. N 

IF ( INOX ( I ) .GT.O)  GO  TO  80 
TEST=A»S(A(1.N*1) ) 

IF(TEST.gT.1.0E-8)GO  TO  99 
30  CONTINUE 

DO  70  1=1. N 

IF  ( INDX 1 I ) .tQ.O)  GO  TO  70 
XIINOXdl  I =A  ( I ,N»  1 ) 

70  CONTINUE 
RtTUHN 

99  MR ITE(2.100) 
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100  FORMAT  120X1  IrtNU  SOLUTION) 
I CK  = 1 
RETURN 
END 


SUBROUTINE  M A T G ( I , J ) 

COMMON/COM  1/00,0 ( 10,11) 

COMPLEX  00 

G ( I , J) = REAL (00) 

G(I,J*l)=-AlMAu(00) 

6(I»1*j*1)=G(I*J) 

G(I*1*J)=-G(I,J*1) 

RETURN 

END 


FUNCTION  SEAST (HH,TT,wW) 

C HH  IS  SIG.  WA/E  HT.  IN  FT,  TT  IS  PERIOD  IN  StC,  WW  IS  FREQUtNCV  IN 

C RAD/SEC.  OUTPUT  SPECTRUM  HAS  UMTS  FT  **2  / (KAD/StC)  . 

COMMON/SSOM/AOO  (80)  ,A10(80>.A01(80),A20(aO),All(aO).A02(dO> 
DIMENSION  F (2) 

H=HM». 3048-4. 016 
T=TT-9.I59 
W=W»»TT/fe. 283185 
IF  (*.GT.0.05)  GO  TO  2 
SEAST  = 0 . 

RETURN 

2 IF  (R.LE.4.0)  GO  TO  3 
SEAST  =0  . 

RETURN 

3 CONTINUE 
N=INT (W/.OS) 

DO  1 1=1,2 
M*N  ♦ I - 1 

1 F ( I ) =A00 (M) *A10 (M) »H* A01  (M) »T*A20 (M) *H*H*A 1 1 (M)  3H*T ♦ A 02 ( M > »T »T 

S = F ( 1 ) ♦ <F (2) -F  ( 1 ) ) » ( w-N*.Ob) “20. 

SEAST=S“HH*“2*TT/101 . 1593 

HETUHN 

END 
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block  data  stAsiow 

COMMON/SSGM/AOO  130  ,A10 (80) , A01  (80)  »A20  I BO  ) ,A1 1 (bO) ,A02(  BO) 

DATA  A00/0.»0.».uo001,.00018,.00133,.00324,.00709,.0l32b..02bl8. 

I .05336,  . I 164  1 , .2603.  .4943,  .8  Job  A,  1 .23195. 1 .59871  . 1 . 7 99S5.  1 . 76253, 
21.56762,1.30231,  1 . 0 7908  . .9  l 784  , . 7 7 733  . .668 1 6,  . 5 7326  . .43533. 

3. 38482..  33 183.. 2828 7..25230. .23205.,? 1658. .2037, . 19481  , . 18371  . 

4. 17350.. 16129. .14752. .14327. .13558. .12091.. 10697.. 09764. .09052, 

5. 08372..  0764,.. 06884,. Ob932». Ob 1 So , .04  350 ♦. 036o0 , . 03037.. 02363. 

6.01831 . . 01466. .01 1 1 7. .00829. .00,bl , ,0039b,. 00283. .0022b,. 00143, 

7. 00057..  00006, -.00041  ,- .00032 ,-. 000 1 2 000  Ob .0 0038 00059 , 

8- . 00077,-. 0009 7,-. 00080. -.00047,-. 00032,-. 00022,-. 000 14, -.00008, 

9- . 00003/ 

OATA  A10/0. 0,0.0, -.00001, -.00004, -.00043, -.00134,  .00255,  .00387. 
1-. 00543, -.004 75, -.000 17,. 0090 1 ,. 02629 ,. 04993 , .06652 ,. ObOOO . .0 3906 , 
2.00467. -.03727, -.06926, -.07963, -.06424, -.05265, -.043j2»-.0326l, 

3-. 01 857, -.01 263. -.009 11, -.00801 ,-. 00336 ,. 00342, . 00539 ,. 00458 ,. 004 , 

4.00652. . 00907. .00923.. 01 084.. 0161 3. .01451 ..01 0,3, . 00d39 ,. 00592, 

5. 00532. . 00714. .00877. .01007. .010 7 7. .01001. .00923. .00750. .00467, 
6.00175,  .00034, .00066,  .00106, .00095, .00090, .00102, .00091  , .00068, 

7. 00036..  00 050.. 000 7 7.. 00093.. 000 73.. 00 027.. 000 18.. 00015.. 000 03, 

8-. 000 09. -.0001 3. -.0001 3, -.000 09, -.000 06. -.00003, -.00001, -.00001, 

90.0. 0.0, 

DATA  AO  1/0. 0,0.0, .00001,. 00003, 0.0. -.00067, -. 0024 ,- .00558 ,- . 00822, 
1-. 01 065,-. 01 169,-. 01241 ,-.00664, . 0 1 278 , . 039 74 , . 06999 , . UH 1 7 7 , . 0558 , 

2.01841 .. 0027, -.00276, -.01 522. -.03524,-. 03485,-. 03 189, -.03983, 

3- . 03554, -.03005, -.02822, -.02864, -.02787, -.02231, -.01716, -.01219* 

4- . 01098, -.01213, -.01061, -.01317, -.02021, -.00812,. 00344, .00783, 

5.01083. . 01 190. .01 1 13. .01021 . .00988. .00930. .01 11,,.01152,.011b4, 

6.01 193. . 01243. .01 189. .01054. .00913. .00785.. 006 74. .00554. .00475, 

7.00422. . 00403. .00345. .00256. .00184. .00129. .00124. .00 120.. 00 109, 

8.00093. . 00073. .00051 . .00027. .00021 . .00023. .00021 ..00 01 7.. 00013, 

9.00007. . 00004, 

DATA  A20/0.,0.,0.,0.,0.,. 00005, .00009, .00016. .00035, .00033,. 00022, 

1.00079..  001 72.. 0041 7.. 00481.. 001 19 ,-. 0066 ,-. 00935 ,-. 00604 ,-. 00044 , 

2. 00 188..  00049.. 00021, -.00021. -.0003, -.00107, -.001 37. -.00 081, 

3.00131 .. 00251 ..00 183.. 00020 ,-.00063,-. 000 76, -.000 8 7 ,-.0006, -.0005, 

4.00013. . 00108. .0005, -.00 001 , .00023, .00042, .00046, . 00045 .. 00052 , 

5. 0003..  00017, -.00002, -.00015, -.000 11*. 00002,. 000 11,. 00008, -.00012, 
6-. 00028, -.000 12,. 000 11,. 00036,. 0004,. 00035,. 0001 3,. 00027.. 00054. 

7. 00058..  00040,  -.0  00 03, -.0  00 14, -.0  001  ,.00001,. 00014,. 0002,. 000 22, 

8. 000 15..  00004, -.00001, -.00003, -.00003, -.00002, -.1)0001/ 

DATA  m11/0.,0.,0.,. 00002, .0002, .00041, .00077, .001 12, .00146, ,00103, 

1- . 00103,— .00667,— .01387,— .02494,— .02849, -.01366, .01256. .024l<», 
2.0251  3,.  01  785..  01  365*.  01  369,.  01  28  7,.  0 1 19,.00914,.0C60<,,.00  4<,1, 

3. 00222,-. 00 30 3. -.00754,-. 00807. -.00403, -.00067,. 00046. .00026, 

4- . 00 086, -.00 096. -.0031 ,-.00 798, -.00544, -.00238, -.00232. -.00222, 

5- . 00222. -.00281  00349 ,-. 00324 ,-. 00292 ,-. 0 0 1 99 ♦-.  ) 0 1 <*6  ,-. 0 00 l 1 , 

6- . 00094 ,-.0004  7. -.0000  7,  . 00031,. 00056,. 000 19,-. OOOcb. -.0  00 7 3, 

7 - .000 75. -.00059,-. 000 10. -.00016,-. 00055.-. U0054,-. 000 3, .U0U3, 

8. 0004 8. . 000 36.. 0. 1016, -.00002, -.0001, -. 0001 l,-. 00 005, .000 06.. 00 01J. 

9. 00015..  0001 4,. 00009,. CO 005/ 

DATA  A02/0.,0.,0.,-. 00004, -.00021, -.000 16.. 00027,. 00 14.. 00 193, 

1 .00 188. . 00082.. 00 042, -.00 032. .00428,. 00436. -.00658, -.02142,-. 01 77. 

2- . 01  loo. -.00-,  11  *.000 16,. 00259.. 0 0845,. 0 081 8,. 00924, .01 304..0 1044, 

3.00776. . 008 19., GO 99b, .00943, .00585, .00222, -.0001 1 ,-.001 39. -.00171. 

4- .00314,-. 00183, .00268, -.00207, -.00614, -.0063, -.00615, -.00599, 

5- .005J2, -.00432, -.0036,-. 0029, -,U03*-,00267,-. 0021 1 ,-.00171 . 


I 


b-. 00145,-. 00095,-. 00043,-. uoooi..ou os i..ouo94, .001 2b. .00134* 

7. 00  1 33. . 0 0 lid.. 00 1 Ob. . 00104, .U004..0007a,.ou0bt'».0005.. 00044, 

’I  8. 000 43, .0003b». 00029, .000 25 ,.000c2, ,0001b.. 00014, .OOOOB..OOOU4. 

"i  9.00002, .00001/ 

END 


i 


I.  < 

I I 

| I 


i 

i 

I 


SUBROUTINE  ZEHU ( A.B22.B44.Bbb) 

COMMON/Nt  *6/F  I AV  . YA  »A0  , SHAY  AV 

C0MM()N/<NE<*b/O«M  (o)»5(6>,SA(6),Y(b),Z(6).NFS.B<b)  ,C05S  lb)  , SING  <b> 
CGMPLE A G22 , 8 44  » t>bb 
00  1 1=1,  NFS 
T = .B488»Ri*5  ( I ) 

P=GAM  (11/57.3 

IF ( I . Lb • 3 ) GO  TO  94 

IF  ( I .(ir  .4)GO  TO  95 

YH  = Y ( I ) -.5*8 ( 1 ) *COS5  l I ) 

ZH=2<I)..5*B<I)»AbS(SING<1)> 

YHR=Yk*CUSS ( I 1 »/H*SlNGl I) 

BB  = H ( i ) 

GOTO  9b 
95  CONTINUE 
88=2. *8  ( I I 
YHH  = -Z ( I ) - B(I) 

9b  CONTINUE 

ABM  3=  ((YRR»8B>**4-YH**4)/(it.*fc!cj) 

B44=944»l . 1 7*T»AMM3*F IAV 
GO  TO  3 
94  CONTINUE 

5 1 = T AN ( H ) 

52  = -Y  l I ) / L (I) 

ALF  = A8S ( (S2-SI ) / ( 1 . *SI*S2) ) 

ALF  = AT  AN  ( ALF ) 

AHM  = SURT  < Y ( I ) »*2-*2 ( I ) **2 ) 

B44=B44*T*AHM**3*CNS(ALF)  *FIAV 
3 CONTINUE 
ALF=ABS (P) 

T=T*CNS(AlF) 

B22=822*T*SrfAYAV 

Bbb  = B6b»T*YA*AVM  ABS  (SX  ( I ) ) ) **3 

i continue 
return 

ENO 


FUNCTION  CN5  l ALF  I 
A=57.3**LF 
IFIA.LT. 40)  UO  TO  1 
CN5  = 1.1 7*S IN ( ALF  ) 
RETURN 

1 CN5=.04b7*A*S  [NlMLF  ) 

RETURN 
ENO 


FUNCTION  S0AMP(W,6> 

DIMENSION  F ( 161 

DATA  F/0. 0 ». 024.. 0*B«. 29B.. 5 7*.. VOS« 1. 1 24, 1. 2 J8.1. 206. 1. 167 » 1.0 71 
•♦.<JBl».B93..B21..74  7,.6b6/ 

T»6*B*d/32.2 
IF(T.GT.O.O)  GO  TO  I 
SUAMP=0. 0 
RETURN 

1 P»T/0. 1*1.0 
N»INT(P) 

IF (N.LT • 151  GO  TO  2 
N*15 

2 C=F (N) * (P-N) * (F (N*l)-F (N) ) 

IF (C.GE.O.O)GO  TO  3 
C = 0.0 

3 CONTINUE 
SDAMP*C**»*B 
RETURN 
END 


FUNCTION  CLALF (A0.A.2PI) 

AOP I = AO/3. 1 4 1 593 

CLALF  = A0*A/ ( AOP I*ZP 1 *SOHT { A*®2*A0PI*tt2) ) 

RETURN 

END 


FUNCTION  BIPLtrt) 

BIPL  = I1.0-.66*H)/(1. 055*3. 7*H) 

IF  (BIPL  .G£.  0.0)  GO  TO  1 

BIPL  = 0.0 

RETURN 

ENO 


SUBROUTINE  FLAPtP.Tl ,T4»T7,T8.T10.T11) 
P2=P*P 

X1*SQHT (1.-P2) 

X2  = AS  IN  (XII 

T1=-X1*(2.*P 2)/3.*P*X2 
T<*  = -X2*P*X1 

T7*-X2*  ( . 12S*P2) ♦ . 125*P*X1*  <7.*2.»P2> 
TBa-Xl* ( 1 .*2.<*P2> /3.*P*X2 
.T  1 0*X  1 *X2 

T11=X2*(1.-2.*P>* ( 2 . “P 1 *X  1 

RETURN 

ENO 
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I 
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SUBROUTINE  TMc.OouMA.CUCK.  St  I 
COMPLtA  61 .C.C6.C3.G.G6.G3.GI .CK.Sc 
Al  s (0.  0.1.0) 

GI=AI*0 

C=1 .O-UI* (. 165/  ( .0«.5*QI)  *.335/  ( .3*01)  ) 

C6=l.0-.361*UI/< .381*01) 

C3=1.0-.2d3*OI/(.54*Ol) 

G= 1.0-0  I *(.236/ (.058 *01)  * . 5 1 3/ ( . 364*0 1 ) * . 1 7 1 / ( 2 .42*0  I ) ) 
G6=1.0-Gl»(.44d/(.29*QI)*.272/(.72S*U1)*. 193/(3. 0*01)) 
G3=  1 . 0-0 1 * ( .6  79/ ( .556*01) *.227/ (3.2*01) 1 
AC=CABS(C) 

PC=ARGO (C) 

AC6=C AdS ( C6 ) 

PC6=AHGD(C6) 

AC3=C AdS ( C3  1 
PC3=ARGD ( C 3 ) 

AG=C AdS ( G ) 

PG= ARGU (G) 

AG6=C AdS ( G6 ) 

PG6=ARG0 (G6) 

A63=CAdS (G3) 

PG3=AHG0(G3> 

IF  (A  .GT.  6.01  GO  TO  1 
AF  = F36 ( A.AC3. AC6) 

PF  = F36( A.PC3.PC6) 

AG  = F36(A,AGJ.AG6) 

PG  = F361A.PG3.PG6) 

GO  TO  2 

1 CONTINUE 

AF  = FGT6(A.AC3.AC6.AC) 

PF  = FGT6 ( A.PC3.PC6.PC) 

AG  = FGTfa ( A.AG3.AG6.AG) 

PG  = FGTo(A,PG3.PG6.PG> 

2 CONTINUE 

CK  = AF  * ( COS ( PF ) ♦ AI *S IN ( PF ) ) 

SE  = AG* l COS ( PG ) * AI*SIN (PG) ) 

RETURN 

END 


function  F36(a.Y3«Y6) 

F 36  = Y3  ♦ (Yb-Y3) /J.OMA-3.0) 

RETURN 

END 


FUNCTION  FGT6(A*Y3.Y6.YC) 

s s (Y6-Y31/3.0 

AA  = 12 . 0* ( Y6-YC  ♦ 3.0*S> 

0 = -36.0* (6. O^S  ♦ Y6  - YC) 
F6T6  = YC  ♦ AA/A  ♦ 0/A*»2 
RETURN 
END 


FUNCTION  AHGO(Z) 
COMPLEX  Z 
XsREAL ( Z ) 
Y^AIMAulZ) 
ARGD=ATAN2(Y.A> 
RETURN 


SUBROUTINE  HULL  I 

COMMON/GH/NUT  . HON  , C A Y , AMC  * OFC  . 1 A ( 2 b , 8 ) ,2A  (2b,H> 

COMMON  Pi  , HP  I » UP  i i TP  I »MU  , MOUt  ,UPrt,CH,HAT  *bUHtUtO«  XT  tUH  f ttidHibUtN 
10E,PDM,V0t.»0Et.,UN,CMEGA,CP.WVh.  lU.DOG,  IG,  XX  (25.  7)  » Y Y ( 2b  » 7 ) »OtL  ( 25  • 
£7  ) ,SNE(25»7)  ,CSE(25»7)  ,FH(7)  , nLOo ( 2b , 7 , 7 ) , YLUG < 2b . 7 . 7 1 .CON  (14,1)  ,C 
JT  (14,14  / »PSIL  (7.7)  .PSI2I7.7)  ,PRa<  7)  ,PRV  (7) 

COMMON/ NE  W/AA  ( 25)  , OX  A (25)  , XCG  . tL  . NST  , HCG.C4  4H.0  ISP  , HHO 
C0MM0N/NEA6/KIAV,  YAWAV.SWAYAY 
hPI=.5*PI 
UP  I = .5*HP 1 
67  MEAD ( 5 . 1 3) GMI N 

WHITE(6,206)GMlN 

HEAD(5»13)FIAV»YAWAV»SWAYAV 

WHITE (6.207) F I A V , Y AWA V , SW AY AV 

YAWAV=YAW4V/b7.3 

FIAV  = FIAV/57.3 

HEAD ( 5 ,20 1 ) NST 

DO  1 I ST= 1 »NST 

HEAD  (b.44) XA  ( 1ST) 

READ (b, 13) (YA(1ST,j),J=1,8) 

Ht AO (5,13) <ZA(ISr,j) 

WRITE (6.205) XA (1ST) 

WHITE (6,36) 

WHITE (6. 13) (YA(IST,J),J=1,8) 

WHITE(6,37I 

1 wRITE(6. 13)  (ZAIIST.J) ,0=1,8) 

00  45  1=1. NST 
45  XA(I)=XA(I)»EL/ZO. 

0AA(1)=.5*XA(2) 

NP=NST - 1 
DO  65  1=2. NP 

65  DXA(I)=.5*(XA(1*1)-XA(1-1) ) 

OX A (NST ) =£L-.S* (XA (NST) *XA (NST-1 ) ) 

00  66  1=1, NST 

66  XA(I)*XCG-AA(I) 

N0N  = 7 

NUT  >8 

00  424  1=1, NST 
DO  424  J= 1 »NUT 

424  ZA(I,J)=ZA<I,j)-ZA(I,NUT) 

NOE»2»NON 
C44M«0.0 
DO  90  I ST  = l ,NST 
C4  = 0. 

00  89  1=1, NON 

XINT  = YA(lST,I»i)-YA  (IST.l) 

Y1NT=2A( 1ST ,1.1 ) -ZA (IST.l ) 

OEL (IST.l) =SUHT (XINT*XINT*VlNf*YlNT) 

SNE (IST,I)=YINT/OEL(IST»I) 

CSE(IST»I)=XINT/OEl(IST,I) 

XX«IST,I)  = .5*(YA(IST,I,1) *YA(IST,  1) > 
YY(IST,I)=a5*(ZA(IST,I»l).ZA(IST,I)) 

89  C4 *X X ( I ST , I ) * ( X IN T * X X ( I S T , 1 ) ♦ Y 1 N T 4 ( Y Y { 1 S I , I ) -hCG)  ) *C4 

90  C44h  = C44M,64.4>»C-,oUXA  ( 1ST) 

GMCALC=C44H*HHU/(2240.*DISP) 

WHITE (6.208IGMCALC 
IF (GMIN.LE.O.O) GO  TO  230 


1 
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CAAH=22AO.°OI5R°IjMIn/RHO 
230  CONTINUE 
SG=- 1 .0 
MO  = 2 
CR  = 0 

DO  300  IST=i.NsT 
DO  301  J= 1 ,NUT 
YA«IST,J)=YA(IbT,J)/EL 

301  ZA ( 1ST, J) =2A( 1ST. J) /EL 
DO  302  J* 1 , NON 

DEL  ( 1ST  , J)  =l)EL  ( 1ST  , J)  /EL 
A X ( I ST , J ) =XX ( 1ST, J) /EL 

302  YY(IST,J)=YY(IST,j)/£L 
CALL  FINO 

300  CONTINUE 

RETURN 

13  FORMAT (8F10. A) 

3b  FORMAT ( lHO.SX.VHASSCISSAS) 

37  FORMAT  (1H0 . 5X  , <iHOROlN  ATES  ) 

AA  FORMAT (F10. 3) 

201  FORMAT ( 12) 

205  FORMAT ( 1 HO . 7HS T AT  ION ,F6 . 2) 

206  FORMAT l/5X»GMIN=»F8. A) 

207  FORMAT ( /5X«F I A V = »F8 . A ,5X« Y Aw AV  = *F 8 . A ,5X*SwA Y AV=*F8 . A ) 

208  FORMAT  < /5X*GMCALC  = *>F8  . A) 

END 
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SUbWOUTINt  HULLO 

COMPLEX  B2  2 » B21.  , b2b  , bA  A . b Ab  * bb2  tdb1*  • B66 
COMMON/GH/NUT  .NUN.CAY , AMC.DFC, YA(2S,8) .ZA (2b. 8) 

COMMON  PI  ,HPI .DPI . TPI .MD.MUDE  .DPH.CR.RA I . SUM , DEG . 1 S T . OH T .hdM.SG.N 
lOE,PI)M,VGL»OE»»DN,UMEGA,CP.*Vh,ID»0OG, IO.XX<25.7).YY(25»7) .DEL (25, 
27) .SNE (25, 7) »CSt ( 25.7) ,FH ( 7) ,cLOo <25.7. 7) . YLOG ( 25 .7,7) .CON U 4,1) .C 
37 ( l A, 14) ,PSIl ( 7.7) .PS12(7, 7) ,PRA( 7) ,PRV (7) 

COMMON/ NEV./X a (25)  .OXA  (25)  . XCG ,EL , NST » HOG .C44H ,D ISP ,RHO 
C0MM0N/NEx2/M.U*A22»822.A2A.d2A.A2b.B26.A‘»L.bH<..AAb.b<.br 
lA62.Bb2.A6A.boA.Abb.b66.EF ( 10) 

C0MM0N/iVtl*3/A22F.  A2AF.A26F,  AAAF.AAbF,  A66F 

0 IMENb 1 ON  EM (3. 25)  .FI  (3,25)  .Hh(3.25) .Hi (3,25) .AM (3. 25) .OF (3,25) 

Q=W**/32.2 

wL=TPI*32.2/W**2 

00  = EL 

CAY=Q-0D 

UN=CAY 

OMEGA=SQRT (UN) 

DO  100  IbT=l»NbT 
DO  100  M0DE=2.3 
GO  TO  (303, 30J.30A) MODE 
303  DO  305  U=  1 , NON 

305  FR ( J) =-SNE ( 1ST, J) 

GO  TO  au 

30 A DO  30b  J = 1 » NUN 

306  FH ( J)  = ( YY  ( IST.u) -hCG/DD) *SNE ( IST.U) *XX ( 1ST, J) *CSE ( 1ST, J) 

HO  CALL  FREQ 

£W<M0DE.I5T)=0. 

FI (MODE. 1ST ) =0. 

HR(MOUE,IST)=0. 

HI (MODE. 1ST) =0. 

DO  Al  1=1, NON 

02  = EXP (CAY»YY ( 1ST,  I ) ) «DEL ( I ST , I ) *32 . 2*DD 
Q3  = CAY *XX ( 1ST  , I ) 

QA=SIN (03) 

05=C0b(U3) 

Q6=SNE ( 1ST , I ) »u5-CSE ( 1ST , I ) «UA 
FI (MODE. 1ST) =F I (MODE. 1ST) -FH ( I ) *u2*U4 
HR  (MODE.  1ST  ) =HM  (MOCE.  1ST)  ♦02<kPR V ( I ) *06 
Al  HI (MODE, 1ST) =H I (MODE, 1ST) .02*PRV( I ) *06 
GO  TO  (50,50,51) MODE 

50  AM (MOUE • 1ST ) =AMC*D0*OD 
OF (MODE. 1ST ) =w*DFC*00*DD 
GO  TO  100 

51  AM ( 1 ♦ I S T ) = 0 • 

DF ( 1 , 1ST ) =0  . 

ER(3.IST)=00«EH(3,IST) 

FI (3. 1ST) =DO»Fl  (3, 1ST) 

HR (3, 1ST) =DU*nR (3. 1ST ) 

HI (3, 1ST ) = U 0 * H I (3,  1ST ) 

DO  52  1=1, NUN 

AM  ( 1 , I ST ) = AM { 1 . I ST ) -SNE ( 1ST . I ) “PKA ( I ) *OEL  ( 1ST  , I ) 

52  OF ( 1 .1ST) =DF ( l , 1ST)  -SNE < 1ST . I) “PHV ( I ) «OEL ( 1ST, II 
AM ( 1 . 1ST) =AM ( 1 , 1ST ) *6 A, A* (00/ A) **2 

DF  ( 1 , 1ST)  =UF  ( 1 , 1ST)  »6A.4*UU«“2/*» 

AM (MODE .1ST) =AMC*D0*»4 
OF  (MOut . 1ST  ) =OFC*uO*»a*i) 


60 


100  CONTINUE 
U*=U/lrf““2 

307  A22=-U»“DF <2, NST) *A22F 

B22=U*AM(2,NST) *822 
A24=-UW“DF ( 1 *NST ) ♦ A24F 
B24=U“AM(1,NST) *624 

A26  = -U»“XA(NST) “OF  (2.NST) *U“U*“ AM ( 2 ,NST ) ♦ A26F 
B26=U*X A ( NST ) “AM ( 2 » NST ) *U“UW»DF (2, NST) *826 
A44=-UW“DF  (3, NST)  *A<»4F 
644=U*AM(3»NST) *844 

A46=-U»“XA (NST) “DF ( 1 ,NST) *U*U**AM ( 1 ,NST) *A46F 
B46=U*XA(NST)“AM< l.NST) *U“UW*DF ( 1 , NST ) *846 
A62=-Uw“XA  ,NST) “DF (2.NST) *A26F 
B62=U*XA (NST) “AM (2. NST) *862 
A64=-U*“XA (NST)  *DF(1,NST)  *A46F 
Bb4=U*XA (NST) *AM( 1 , NST ) *864 

A66=-U'»“XA(NST)  **2“UF  (2.NST)  ♦U*UW*XA  (NST)  *AM(2»NST)  *A66F 
B66=U*XA (NST) *“2“AM (2. NST) *U“U*»XA (NST > “DF (2. NST) *066 
U2=2.“U/W 
EF ( 1 ) =U2“H I (2. NST) 

EF (2) =-'  ?*HH ( 2 • NST ) 

EF (3) =U2*MI (3. NST) 

EF (4) =-U2*H8(3.NST) 

EF (5) =EF (1)“XA(NST) 

EF  < 6 ) =EF (2)*XA(NST) 

308  CONTINUE 

DO  103  IST=1»NST 
XOX=XA(IST)“DXA(IST) 

XDX2=X a ( 1ST) ““2*0XA ( 1ST ) 

02  = 2 .“DXA ( 1ST ) 

A22=A22*AM(2*IST) “DXa(IST) 

822=822 *DF (2,1  ST ) “DXA  ( 1ST) 

A24=A2a*AM (1,IST)“UXA(IST) 

B24=B24,0F(1,IST)“DXA(IST) 

A2b=A26*AM<2, 1ST) “XOX*UW“DF  ( 2 , I ST ) “DX A ( I ST ) 

826=826 *0F (2 , I ST ) “XDX-U“AM ( 2 , I ST ) “DXA (1ST) 

A44=A44*AM(3,IST) “DXA< 1ST) 

B44  = 8«,4*0F(3*IST)“DXA(IST) 

A46=A4b*AM ( 1 , 1ST ) “XDX*U«*DF ( 1 , 1ST ) “DXA ( 1ST ) 

B46=B4b*0F ( 1 , 1ST) «XDX-U“AM( 1 , 1ST) “DXA (1ST) 

A62=A62*AM(2» 1ST) *XDA-UW*DF (2, 1ST) *DXA< 1ST) 

862=862 *0F (2, 1ST ) “XOX*U*AM (2, 1ST ) “DXA ( 1ST ) 
A64=A64*AM(l,IST)*ADX-Uto*DF(l,IST)*DXA(IST) 

864=B64*DF ( 1 • 1ST ) “XDX*U“AM ( 1 , 1ST) “DXA (1ST) 

A66=Aoo*AM (2. 1ST) “XDX2*U*UW*AM (2, 1ST) *DXA (1ST) 

866=8*6 *UF (2, 1ST) “XDX2 *U*UW*DF (2, 1ST) “DXA (1ST) 
EF(l)=tF(l)*D2“(EH(2,IST)*hR(2«IST)) 

EF (2)=EF (2) *02“ (FI (2, 1ST) *H I (2, 1ST) ) 

EF ( 3 ) =tF ( 3 ) *02* (ER ( j, 1ST) *HK(3,IST) ) 

EF (4) =EF (4) *D2*(FI (3, 1ST) *HI (3, 1ST) ) 

EF (5) =tF (5) *02“ (XA (1ST) “(EH (2, 1ST) ♦HR(2,1ST> ) *U“hI (2,IsT)  /«*> 
103  EF (6) =EF (6) *02“ ( A A ( I ST ) “ ( F I ( 2 . I ST ) *H I ( 2 . 1 ST ) ) -U*hrt ( 2 , 1ST ) /w ) 
HETUHN 
END 
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SUBROUTINE  FINU 

COMMON/GH/NUT .NON# CAY .AMC.UFC.YA (25.B) ,Za < 2b.nl 

COMMON  PI ,MPI .UPI .TPI.MU.MOOE.UPH.tR.KAT.SuK.UEG. Isr.uRT.HhM.SG.N 
IOE«PDM.VOL»OE*.UN,OMEGA.CP."VM. ID.UOG. IG. AX  <2S.  71 .YY<2b«  71  ,utL (2b. 
27 1 »SNt  (25*7)  »CSE  ( 2S  « 7 1 *FH ( 7 1 . BLOG  ( 25. 7 , 7 ) . YLOG  ( 2b . 7 . 7 ) . CON  ( 1 A » 1)  ,C 
3T  <14.141  ,PSIl  (7*71  .PS  12 (7*7)  .PRA(7)  .PKV<  7) 

00  1 1=1. NON 
XMl=AX<IST.Il-YA<IST.l) 

YM1  = YY ( 1ST . I ) A ( 1ST  . 1 ) 

XPI=Aa<IST,II*YA<IST,1) 

YPl»YY<IST.I)*ZACIST.l) 

FPR1=.5*ALOO< AM1**2»YM1**2> 

FPL1=.5*AL0G<AP1*»2*YM1**2) 

FCR1=.S»AL0G(XM1**2*YP1**2) 

FCL1=.5*ALOO(XP1**2-YP1«*2I 
APR1=aTAN2(YM1 ,XM11 
APL 1 =AT  AN2 ( YM1 .XP1) 

AC«1=ATAN2 ( YP1 ,XM1 ) 

ACLI=ATAN2(YP1,XP1) 

00  1 J= 1 » NON 

XM2  = XA(IST.I)-YA(IST»J-»1) 

YM2=YY( IST,I)-ZA(IST,J*1I 
XP2  = XX ( 1ST  » I ) ♦YA(IST.J*1) 

YP2  = Y Y ( 1ST . I ) *ZA ( 1ST » J*  1 1 
FPH2=.5*AL0G ( AM2**2*YM2**2) 

FPL2=.5*ALCG(AP2»«2*YM2**2) 

FCM2*.5*AL0G(XM2**2*YP2**2> 

FCL2=.5*AL0G( AP2**2*YP2**2» 

APH2=ATAN2( YM2.XM2) 

J1=J*1 

IF<XM2.GT.0.UIGO  TO  4 
IF(Jl.GT.I)  GO  TO  b 
IF (YM2.LT.0.0) APR2=APH2*TPI 
GO  TO  5 

IF { YM2.GE.0.0) APR2=APR2-TPI 
IF (YP2.LT.0.0)  GO  TO  4 
ACR2=-PI 
GO  TO  3 
4 CONTINUE 

ACR2=ATAN2(YP2.XM2) 

3 CONTINUE 

ACL2=ATAN2 (YP2.AH2) 

APLZ=ATAN2(YN2«XH2) 

SIMJ=SNE(IST.I)*CSE(IST.J)-SNE(1ST.J)*CSE(IST,I) 
CIMJ=CSE<IST.II*CSE<IST,J) *SNt ( 1ST . I ) *SNE < 1ST . Jl 
SIPJ=SNE(IST.l)*CSE(IST.J) *SNE < 1ST , J I *CS£ < 1ST . X I 
CIPJ=CS£ (IST*I)“CSc(IST*J) -SNc ( 1ST . 1 ) *SNt ( 1ST. Jl 
OPNR=SIMJ* (FPM1-FPh2) .CIMj*  UPk1-APR2) 

RR  PPHaCSE  t 1ST. J) * (XM1»FPRI-YM1 ® aPM 1 - AM l-XM2*F PN2* YM2* APH2* 

1XM2) *SNE ( 1ST. Jl *{ YM1*FPN1»AM1*aPM1-YM1-YM2*FPN2-AM2»APH2*YM2) 
0PNL=SIPj*(FPl2-FPL1) *CIPJ» (APL2-APL1) 

PPL=CSt< IST.JI *<XP2*FPL2-YM2*APL2-XP2-XP1*FPL1*YM1*APL1 » 

1 XP l » ♦ SNt ( 1ST. J) * ( YMi*FPLl»XPl*APLi*YM2-YM2*FPL2-AP2*APL2-YNl I 
DCmR=SlPj*  (FCh1-*-C«2)  *UPo*  ( ACM1-ACR2) 

PCP=CSE  ( 1ST.  J|  • ( AMl*FCRl-YPl*ACrU-XMl-AM2*FCK2*YP2*AC:-i2* 

1XM2) ♦SNE ( 1ST, j) * ( YP2*FLH2*XM2*ACM2.YP1-YP1 »FCK 1 -AM  1 »AC* 1 - YP2 ) 
DCNL=SIMJ* (FCL2-FCL1 J *C IM J* ( ACL2-ACL 1 I 


% 


PCL*CSE  ( 1ST  *U) * (XP2*FCL2-YP<;*ACL2-AP2-AP1  »FCL  1 OH  1* ACL  1 *XP 
11) *SNE  ( IST*J)  »<YP2*FCL2*AP2*ACl2-YP2-YP1*FCL1-XP1«ACH  *YP)  ) 
BLOG ( 1ST, 1*0)  *OPNR*SG*DPNL-UCN*<-Sli»OCNL 
YLOG  ( 1ST  * I . J ) =PPR*SG*PPL-PCH-SGi*PCL 
IF (d-NON) 2*1 *1 
2 XM 1 = XM2 

YH1*YM2 
XP 1 *XP2 
YP 1 *YP2 
FPR 1 *FPP2 
FPL 1=FPL2 
FCR1=FCR2 
FCL 1*F  CL2 
APR  1 »APR2 
APL 1= APl2 
ACR1»ACR2 
ACL1*ACL2 
l CONTINUE 

RETURN 
ENU 
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SUBROUTINE  FHtJ 

COHWON  HI  ,hHI  ,UR1  . TH1 .MU,MOOE.uHh.CH,HAT,SUR.IItG»IST.OwT,HdM,So»N 
lOE.PDM,vOL.Lt*»uN,OMEGA,CP»*Vh,Il>,OOG»lG.XX  (2S.7)  ,YY  (2S, 7)  ,LcL(25 
2,7)  .SNt  (25.  7)  .LSE  (2b. 7)  ,FR(7>  .BLOG  ( 2S  » 7 » 7 ) . YLOG 1 2S , 7 » 7 ) .CONI 
j»CT<l*,l*)  ,Fiil  (7.7)  ,PSI2(7,7>  ,HHa (7)  ,HHV  ( 7) 

COMMUN/GH/NUT  ..mON.CAY  .AMC.UFC.YA  (2S.9) ,2A (2S.B) 

1U  UO  1 1=1. NON 
NI=NUN» I 

con<i»i)=o. 

CON  (N  I . 1 ) =OMfcGA'ir«(  I) 

XRl=UN» (XX (IST.i)-YA(IST.l)) 

YH1=-UN»( YY ( 1ST . I) *2A ( 1ST . 1 ) ) 

XL1=UNMXX (1ST. I) ♦ Y A ( I ST  » 1 ) ) 

YL 1 =YH 1 

CALL  LAV i 0 (XRI.YH1.EJ1.CXR1.SXR1.HAR1.RBR1 » Ck  1 .SRI ) 

CALL  LAVIO  (XL! . YL 1 » £ J 1 tCXLl .SXLi , HAL  1 .HdL i.CLl.SLl) 

00  1 J=  1 » NON  / 

NJ=N0N* J 

XR2=UN* ( XX (IST,I)-YA(IST»J*1)) 

YR2=-uN* ( Y Y ( I ST  ♦ I ) *ZA(IST,J*1) ) 

XL2=UN* (XX ( 1ST, 1) *YA( 1ST, J*l> ) 

YL2=YH2 

CALL  LAVIO  (XH2.YR2.EJ2»CAH2,SXR2.HAR2,RBR2»tH2»SH2) 

CALL  ,'AVIU  (XL2.Y!_2,tj2,CXL2»  SXL2 .HAL2.RBL2  .CL2.SL2) 

SIPJ=SNE< 1ST, I) »CS£ ( 1ST .J) + SN£ ( 1ST , J) »CSE ( 1ST , I > 

C IPJ=CSE  ( 1ST,  I ) »CSE(IST,J)-SNt(IST.I)*SN£(IST,vJ) 

SIMJ=SN£  < 1ST, 1 1 *CSt ( 1ST. J) -SNt ( laT, J) acS£ ( 1ST, I ) 

C IM J=CSE ( 1ST , I )>CSE l 1ST , J) .SNt ( IST,I)*SN£ ( IST.J) 

CT ( I. J) =bLOG( (ST . I . j ) *2.o (SIHu<* (CH i -CR2 ) -C IP J* ( Sk 1 -SK2 ) -SG* ( S I 
IMJMCLI-CL2)  -C i MJ*  ( SL  1 -SL2 ) ) ) 

PS  II ( I , J) =YLOG (IST,I,J)*2./LN« (SNE (1ST, JI*(HAR1-KAR2)+CSE (IST.J 
1 )* (R8R1-RBR2) »SG<MSNE ( I ST . J) * (HAL  1 -RAL2) *CS£ ( 1ST . Jl * (RBL2-HBL1 I ) ) . 
CT(M.NJ)=CT(I,J) 

CT ( I.NJ>=TPI<*(t J2»(SXR2*CIPJ-CXH2*SIPj)-EJl*(SXHi*CIPJ-CX 
IR1*SIHJ)-SGM£u2*(SXL2*CIMJ-CaL2*SIMJ>-£J1MSXH*CIMJ-CXl1 
2*SIMJ) ) ) 

PS  1 2 ( I » J) aTP I /uN° ( E J1 * ( SXK1 *CSt ( 1ST  » J) -CXh l *SNE ( 1ST  » J) > -E J2* 

1 (SXH2*CSE( IST.J )-CxH2*SNE (IST.J) ) -SG* ( £ J 1 * ( SAL  1 *CSE ( I ST . J ) ♦CXLl'SN 
2E(IST»J))-Ej2*(SXl2*CSE(IST,J)»CXL2*SNE(IST,J) ) ) ) 

CT(NI»J)=-CT(I ,NJ) 

IF ( J-NON) 7,1.1 
7 XR1=XM2 

YR1=YR2 
XL  1 “XL2 
YLI=YL2 
E J 1 =E J2 
CH1=CH2 
SR1=SH2 
CL1=Cl2 
SL  1 =SL2 
HAH l =H AH2 
RDH 1 =HbH2 
HAL1=MAL2 
HBL 1 =HBL2 
CXHI=lXH2 
SXR l =SXH2 
CXL 1 =CXL2 
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SXL  1 = S>  AL2 

1 CONTINUE 

CALL  MAT  INV  (CT  .NOE.CON.  l .DOG.  ICI 
GO  TO  12.6)  *10 

2 00  3 1=1. NON 
PHA<n=0. 

PHV ( I ) =0 . 

DO  A J= 1 . NON 
NJ=NON* J 

PNA<I)=PPA<i> ♦CONIJ,1)*PS12(I»J)-CON(NJ.1>*PS11 U.J» 

4 P«V  (1) = PH  V ( I ) *CON( J» 1 ) *PSI  1 II  * J) *CON(NJ. 1 > “PS12 1 1 » J1 
PHA ( I ) =OMEGA*PHA( l) 

3 PRV ( I ) =QMEGA®PRV  < 1 ) 

AMC=0. 0 

DFC  = 0 .0 
00  5 1=1. NON 

AMC=AMC*PRA ( I I *DEL (IST»I1*FHCI> 

5 OP  C=DFC  »Prt V (I ) <*DEL ( 1ST • I ) *KR ( I ) 

AMC=2 • 0®AMC 

OFC=2.0*DFC 

AMC=AMC/ON 

OFC=DFC/UN 

6 RETURN 
END 
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OAVI  - COMHO T A i ION  Or  FHtOOENC  Y DEPENDENT  PakTS  OF 

2-0  POTENTIALS  ANO  KERNELS 

SUBROOT INE  OAVID! A,Y,t.C*S.RA,MbfCIN,aON) 

AT  = A T AN2 ( X * Y ) 

AHG  = A T-l .5707S63 
E=EXP(-Y) 

C = COS  (A) 

S = SIN  (A) 

R=X*»2* Y**2 
TEST=0. 00001 
IF(R.LT.1.0>  GO  TO  5 
TE5T=0.1*TEST 
IF (R.LT .2.0)  GO  TO  5 
TEST=0.1*TEST 
IF (R.LT.4.Q)  GO  TO  5 
TEST  = 0. 1*TEST 
5 AL=O.S*ALOG (H> 

SUMC=0.S7721bbo»AL*Y 
SOMS=AT*X 
TC  = Y 
TS3  A 

DO  1 K= 1 » 50  0 
T0=  TC 
CO  A =K 
CAY  =K ♦ 1 

FACT  = COA/ (CAY»Cay ) 

TC=FACT*(Y*TC-a»TS) 

TS=FACT* ( Y*TS*A»TO) 

SOMC=SOMC*TC 
SUMS=SOMS*TS 
IF(K.GE.SOO)  GO  TO  3 

IF( (ABS1TCI ♦Aba(TS) I .GT.TEST)  GO  TO  1 

3 CIN=EMC«SUMC*S*SOMS) 

SON  = E*  <S*SUMC-C»SUMa) 

R A* AL-C  IN 
RS*ARG«SON 

GO  TO  4 
1 CONTINOE 

4 HETOHN 

END 


_ 


66 


ooo  non  o o o ooooonnooooooooono 


SUBPOUTImF  MATINV  < A.NR.B.NC.OFTFRM, ID) 

MATRIX  INVERSION  WITH  ACCOMPANYING  SOLUTION  Of  LINEAR  EQUATIONS 
PIVOT  METHOO 

FORTRAN  IV  SINGLE  PRECISION  WITH  ADJUSTABLE  OlMENSION 
FEBRUARY  1B66  S GOOD  OAVID  TAYLOR  MODEL  BASIN  AM  MATA 
WHERE  CALLING  PROGRAM  MUST  INCLUDE 

DIMENSION  A (NR*Nrt) * fl(NR*NC>»  INDFX(NR,3> 

N IS  THE  ORDER  OF  A 

M IS  THE  NUMBER  OF  COLUMN  VECTORS  IN  B (MA Y BE  0) 

DETERM  WILL  CONTAIN  DETERMINANT  ON  EXIT 
ID  WILL  BE  SET  BY  ROUTINE  TO  2 IF  MATRIX  A IS  SINGULAR 
1 IF  INVERSION  WAS  SUCCESSFUL 
A THE  INPUT  MATRIX  WILL  BE  REPLACED  BY  A INVERSEE 

B THE  COLUMN  VECTORS  WILL  BE  REPLACED  BY  CORRESPONDING 

SOLUTION  VECTORS 
INDEX  WORKING  STORAGE  ARRAY 

IF  IT  IS  DESIRED  TO  SCALE  THE  DETERMINANT  CARD  MAY  BE 

OELETEO  AND  DETERM  PRESET  BEFORE  ENTERING  THE  ROUTINE 

EQUIVALENCE  (IPOW.JROW),  ( ICOLUM, JCOLUM) * < AMAX * T*  SWAP) 
OlMENSION  A (NR*NR) ,9<NR*NC> *INDEX(30«3) 

Nl=NP 

M1=NC 


INITIALIZATION 

N=Nl 

M=M1 

DETERM  = 0.0 
DO  20  J=1,N 
20  INDEX ( J»3)  = 0 
DO  550  1=1, N 

SEARCH  FOR  PIVOT  ELEMENT 

AMAX  = 0.0 
00  105  J=1,N 

IF ( INDEX ( J,3) -1 ) 60,  105,  60 
60  DO  100  K=1,N 

IF(INDEX(K,3)-1)  BO,  100,  715 
BO  IF  ( AMAX  -ABS  <A(J»K>)>  85,  100,  100 
85  IROW® J 
ICOLUM  =K 

AMAX  = ABS  ( A ( J,K ) ) 

100  CONTINUE 
105  CONTINUF 

INDEX (ICOLUM, 3)  = INDEX ( ICOLUM, 3)  ♦! 

INDEX (I«l)=IROW 
INOFX ( I , 2) = I COLUM 

INTERCHANGE  ROWS  TO  PUT  PIVOT  ELEMENT  ON  OIAC-ONAL 

IF  (IROW-ICOl.UM)  1 AO , 310,  140 
140  DFTFPMs-DETFRM 
DO  POO  L=1.N 
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SWAP=A < IPOW.L) 

A ( I ROW .L ) =4(1 COLUM.L ) 

200  A ( ICOLUM, L)=SWAP 

IE(M)  310,  310,  210 
210  00  250  L*l*  M 
SWAP=R<  TRON.L) 

R(IPOW.L)=R( ICOLUM ,L ) 

250  R<ICOLUM,L)=SWAP 
C 

C DIVIDE  PIVOT  ROW  BY  PIVOT  ELEMENT 

C 

310  PIVOT  =A  ( ICOl.UM*  ICOLUM) 
DETEPM=OETERM*PI VOT 
330  A ( ICOLUM , I COLUM) =1.0 
00  350  L= 1 ,N 

350  A ( I COLUM.L)  =A  ( I COLUM.L)  /PIVOT 
IE(M)  380,  380.  360 
360  00  370  L=1,M 

370  R ( ICOLUM, L) =R ( ICOLUM ,L ) /P I VOT 
C 

C PEDUCE  NON-PIVOT  ROWS 

C 

380  00  550  L 1 = 1 ,N 

IE(Ll-TCOLUM)  400.  550.  400 
400  T=A(L1. ICOLUM) 

A(L1.ICOLUM>=0.0 
00  450  L=1.N 

450  A(L1,L) =A(L1,L)-A(IC0LUM,L)*T 
IFIM)  550.  550.  460 
460  00  500  L=1.M 

500  B (Ll,L)=fl(Ll.L)-B( ICOLUM, L)*T 
550  CONTINUE 
C 

C INTERCHANGE  COLUMNS 
C 

00  710  1=1, N 
L=N.1-I 

IE  (INOEX (L,1)-INDEX(L,2)  ) 630,  710,  630 
630  JROW= INOEX (L , 1 ) 

JCOLUM=INOEX(L,2) 

00  705  K=1.N 
SWAP=A (K , JROW) 

A (K , JROW) =A (K , JCOLUM) 

A (K  . JCOLtJM)  =SWAP 
70S  CONTINUE 
710  CONTINUE 

00  730  « = l.N 

IE ( INDEX (K . 3)  -1)  715.720,715 
720  CONTINUE 
730  CONTINUE 
10  = 1 
810  RETUPN 
715  10  = 2 

GO  TO  810 
ENO 
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